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SUMMARY 


This research under NASA Grant NGR 03-002-044 was initiated in 
1965, and included the theoretical research for the development of an 
effective methodology for designing specified reliabilities into 
mechanical components, and experimental research to develop three fatigue 
reliability research machines which can apply a reversed bending moment 
combined with a steady torque to round, rotating, ungrooved and grooved 
specimens. 

Phase I of the experimental research program, initiated in 1967, 
included the generation of distributional cycles-to-failure versus 
alternating bending stress (S-N) diagrams and of distributional Goodman 
strength diagrams for specimens made of AISI 4340 steel, R c 35/40 
hardness, and having a circumferential groove which provides a theore- 
tical stress concentration factor, K^, of 1.42. 

Phase II of the experimental research program was identical to 
that of Phase I except that the specimen groove provided a theoretical 
stress concentration factor of 2.34, and was initiated in September 1970. 

Phase II results are compared with Phase I results in this report 
and the effects of superimposing a steady torque on reversed bending on 
the distributional S-N and Goodman diagrams are presented, as well as the 
effect of different K^'s. Such distributional data has to be generated 
to enable the designing of specified, target reliabilities into components. 


xiii 


A methodology and a computer program were developed for the 
generation of finite-life, distributional Goodman diagrams. This 
methodology provides the capability for optimizing a design to achieve 
a target reliability for a specified component life. 

A FORTRAN computer program was developed to estimate the parameters 
of the three-parameter Weibull distribution representing the cycles-to- 
failure data, and to perform Chi-Squared and Kolmogorov-Smimov goodness- 
of-fit tests. The program also calculates the predicted component life 
for a specified reliability. 

Also a study was made of the cumulative fatigue theory found in 
the current literature, and a number of methods for making cumulative 
fatigue reliability predictions were explored. The most promising method 
of conditional probabilities is proposed for further study and verification 
through a cumulative fatigue test program. 


xiv 


1 . INTRODUCTION 


1 . 1 Background 

The research under NASA Grant NGR 03-002-044, initiated at The 
University of Arizona in September 1965, included theoretical research 
for the development of an effective methodology for designing specified, 
target reliabilities into mechanical components and experimental research 
to generate distributional, statistical S-N and Goodman diagrams to 
provide the design data needed in support of this theoretical research. 

During the first reporting period a basic methodology for design 
by reliability in combined-stress fatigue with time dependent strength 
distributions was developed. Mathematical methods in dealing with the 
functions of random variables involved were discussed. Concurrently, a 
supporting experimental fatigue research program was planned. It was 
found that there were no research machines available which could apply to 
round, rotating test specimens a combination of a reversed bending moment 
and a constant torque. As a result three machines, similar in principle to 
the Mabie and Gjesdahl [1] test machines were designed and built at The 
University of Arizona. A complete discussion of the design and development 
of the test machines was given in the report to NASA CR - 72836 [2], 

During the second reporting period, the operational capability of 
the first machine was proven and two additional machines were fabricated. 

The design and development of, and the results obtained from, these three 
machines were presented in the report to NASA CR - 72838 [3]. During the 


1 


2 

same period calibration constants were determined for each machine so that 
the nominal bending stress and the shear stress in the specimen groove can 
be calculated. The calibration procedure, data, analysis, and constants 
were presented in the report to NASA CR - 72839 by Kececioglu and McConnell 
[4], Calibration constants are needed because the bending and shear 
stresses cannot be monitored directly. Strain gages mounted in a specimen 
groove would be destroyed when the specimen failed. Hence it was necessary 
to mount strain gages on the specimen holders instead of on the specimens 
to monitor the bending and shear stresses. This necessitated the 
determination of calibration constants and equations to relate the strain 
at the strain gages to the nominal stresses in the specimen groove. 

Another report by Kececioglu and Smith, NASA CR - 72835 [5], presented 
all of the experimental data generated up to June 30, 1970. Included were 
the reduction of the data, the application of the design by reliability 
methodology, the conversion of the cycles-to-failure data to stress-to- 
failure distributions, and the development of three-dimensional Goodman 
fatigue strength surfaces;, which is the ultimate form of the reduced data 
for direct use by designers. During this period computer programs for use 
in the reduction of the data were developed and refined. The most important 
were program STRESS and program CYTOFR. Program STRESS calculates the 
bending and shear stresses applied to each specimen, the ratio of 
alternating to mean stress, and the mean and standard deviation of each 
stress for each test series. It also calculates the cycles-to-failure from 
the times-to-failure data recorded during the fatigue tests. Program 
CYTOFR calculates the mean and standard deviation of the cycles-to-failure 
for the normal and the lognormal distributions which approximate the true 
distribution of the data. It then calculates the coefficients of skewness 

) 
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and kurtosis and applies the Chi-Squared and the Kolmogorov-Smirnov 
goodness-of-fit tests to determine which distribution provides a better 
fit to the data. 

1.2 Experimental Research 

Schematic diagrams of the NASA Complex-Fatigue Research are shown 
in Figs. 1 and 2. A test specimen is subjected to a bending moment by 
weights hung at the end of a lever arm. Torque is applied through the 
Infinit-Indexer which rotates shaft A with respect to shaft B and holds 
the relative position of the shafts. Two four-arm strain gage bridges are 
mounted on the toolholder in the positions shown in Figs. 3 and 4. The 
output of one bridge is proportional to the strain resulting from the 
alternating bending stress, and the output of the other bridge is 
proportional to the shear strain resulting from the steady torque. The 
numerical relationship between the strains measured at the toolholder and 
the nominal stress in the specimen groove is established through the 
calibration program discussed by Kececioglu and McConnell [4] . 

Two types of tests are conducted: One to determine the cycles-to- 

failure distribution at a given alternating bending stress level and 
bending to shear stress ratio. The other to determine the endurance 
strength, or stress-to-failure, distribution. The steps involved in each 
type of test are summarized in Fig. 5. A bending stress level and a 
stress ratio are selected from the overall test program and 

★ 

assigned to one of the three research machines. A PDP-8 computer program 
is run to determine the corresponding shear stress and the number of 
Visicorder divisions required to represent the bending and shear stresses. 
The Visicorder records the amplified outputs of the strain gage bridges. 

T 


See Appendix A. 
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The test sequence begins with the installation of a specimen in the 

toolholder collets, with the groove centered between the collets. With 

the collet on the strain gage side tightened, the instrumentation is 

zeroed and calibrated. After tightening the other collet, weights are 

added to the bending- load arm, and the torque is applied to obtain the 

desired stresses as indicated by the number of divisions on the Visicorder. 

At this time a microswitch, which stops the machine when the specimen 

fails, is set and an interconnected clock is set to zero. The machine is 

then started, and run at constant speed until the specimen fails. The 

time to failure is recorded and subsequently used to calculate the number 

of cycles to failure for that specimen. After running 35 specimens at 

the specified nominal bending stress level and stress ratio, a data deck 

* 

is prepared for program STRESS to be run on a CDC-6400 Computer. The 
program inputs include the time to failure, Visicorder resistances and 
divisions used for its calibration, and the divisions recorded during each 
specimen for the bending and shear stress levels. The program calculates 
the achieved bending and shear stresses, stress ratio, and cycles to failure 
for each one of 35 specimens used in each test series. Then it calculates 
the statistical mean and standard deviation of the bending stress, shear 
stress, and stress ratio achieved in each test series of 35 specimens. 

The cycles to failure of each specimen becomes the input into program 

★ * 

CYTOFR for analysis of the statistical distribution of the 35 cycles-to- 
failure data. The normal and lognormal distribution parameters are 
determined, as well as .the skewness, kurtosis and the Chi-Squared and 
Kolmogorov-Smirnov goodness-of-fit test results. The lognormal distribution 
parameters are then used to construct the distributional S-N diagrams. 

* See Appendix B. 

** See Appendix C. 
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The Weibull distribution having gained favor recently in fatigue 

studies, it was decided to determine the parameters of the Weibull dis- 

* 

tribution that represents the cycles-to-failure data. PROGRAM WEIBULL 
is now used to accomplish this and to see which one of the three 
distributions (normal, lognormal, and Weibull) represents the cycles- 
to-failure data best. 

The staircase method of testing is used to determine the stress- 
to-failure, or endurance strength, distribution parameters. The 
endurance strength is taken to be normally distributed, and is plotted 
along with the cycles-to-failure distributions to complete the distri- 
butional S-N diagrams. The staircase results are also used to construct 
the very valable distributional Goodman diagrams. 

The experimental research for this reporting period consisted 
of completing the test program planned for Phase I, consisting of cycles- 
to-failure testing at a stress level of 65,000 psi at the stress ratio 
of 0.44 and endurance strength testing, to complete the S-N and Goodman 
diagrams for Phase I research. 

In addition the Phase II research was undertaken. The 
experimental research of Phase II was a continuation of the research 
performed in Phase I, but with new test specimens. 

The Phase I specimens, shown in Fig. 6 were shafts made of 
AISI 4340 steel, MIL-S-5000 B, Condition C-4, Rockwell C 35/40, grooved 
to provide a theoretical stress concentration factor of 1.42. The 
specimens for Phase II, shown in Fig. 7, are identical to the Phase I 
specimens except that they have a different groove geometry to provide 


* See Appendix D. 
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a theoretical stress concentration factor of 2.34. The results of this 
research are presented here and compared. 

Much valuable fatigue reliability design data has thus been 
generated experimentally in support of the theoretical methodology for 
designing specified reliabilities into rotating components subjected 
to combined reversed bending and steady torque. 

On the theoretical side a methodology for generating finite 

life distributional Goodman diagrams was developed. In addition a 
promising method to calculate the reliability of rotating components 
subjected to cumulative fatigue loads was developed and is presented 
here. 


2 . 


EXPERIMENTAL RESEARCH 


2.1 Phase I Research Completion. 

The results of Phase I experimental research accomplished prior to 
this period were reported by Kececioglu and Smith [5] . They included 
endurance tests at stress ratios of <*>, 3.5 and 0.83. Cycles-to-failure 
tests were accomplished at two alternating stress levels for the stress 
ratio of 0.44: 69,000 psi and 60,000 psi. The endurance run for the 

stress ratio of 0.44 was also begun during the previous reporting period 
but was not completed. 

Another cycles-to-failure test series for the stress ratio of 0.44 
with eighteen (18) specimens was run on Machine No. 1 at a nominal bending 
stress level of 65,000 psi. Program STRESS was run on the CDC-6400 
Computer to determine the nominal bending and torque stresses in the groove 
of the specimen and the resulting stress ratio achieved in these tests. 

The results are summarized in Table 1. From Table 1 it can be seen that 
the achieved standard deviations are approximately two percent of the mean 
values of the stresses and the stress ratio. The run was, therefore, under 
control, and the test results were considered acceptable. 

The cycles-to-failure calculated by program STRESS for each specimen 
were used as inputs into program CYTOFR for the analysis of the distribution 
of the data. A summary of the outputs from the program is shown in Table 2. 


7 
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A comparison of the goodness-of-fit results obtained for the 
normal and lognormal distributions shows the following: 

1. The K-S test does not reject either distribution since the 
maximum D value is less than the allowable in each case. 

2. The coefficient of skewness is significantly larger than 
zero for the log cycles than for the straight cycles; where 
the coefficient of skewness of the symmetrical normal dis- 
tribution is zero. 

3. The coefficient of kurtosis for the data fitted to a normal 
distribution is smaller than the value of 3.0 for the normal 
distribution, and the value obtained for the lognormal dis- 
tribution is greater than 3.0. 

Based on the overall statistics and the conformance of previous 
results, the lognormal distribution was chosen. The results are plotted 
in the S-N diagram of Fig. 8. Figs. 9, 10, and 11 give the S-N diagrams 
-for the previous research results of purposes of completeness and 
comparison. Table 3 summarizes the results used to obtain Figs. 8 thru 
11. Table 4 summarizes the results of the normal distribution fit to 
all Phase I cycles-to-failure data and Table 5 of the lognormal dis- 
tribution fit, for purposes of completeness and comparison. 

Tests to determine the endurance strength for the stress ratio 
of 0.44* were run on Machine No. 2. The staircase method [5, p. 19] 
was used with specimens tested to 2.5 x 10^ cycles for success. 

* The stress ratios actually achieved in these tests averaged out to 
r g = 0.45, consequently these endurance strength results are reported 
as being for r g = 0.45. 
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The results for 37 valid data points are shown in Fig. 12. Using the 
equations presented in Mood and Dixon [8, p. 114] the mean and 
standard deviation of the endurance strength distribution were calcu- 
lated, as shown in Table 6, giving an endurance strength mean of 
49,600 psi and a standard deviation of 3,700 psi. The results obtained 
from the cycles-to-failure tests and the endurance tests were used to 
complete the S-N diagram shown in Fig. 8. 

The endurance strength distribution parameters were applied to 
equation [5, p. 37]. 

, 1/2 


S = S (1 + — 
r a v - 2 


-) 


( 1 ) 


and 


1/2 


J S °S ^ + - 2 

r a r 


( 2 ) 


to obtain the parameters of the distribution along r g = 0.45, with 
the following results: 


= 120,900 psi, 

and 

Og = 9,000 psi. 

r 

The incorporation of this r s = 0.45 strength distribution into 
the Goodman diagram shown in Fig. 13 completed the distributional 
Goodman diagram for 2.5 x 10^ cycles of life and the Phase I experi- 
mental research program. Figures 14, 15, and 16 give the staircase 
test results for r g = «, 3.5, and 1.0 respectively; and Tables 7, 8, 
and 9 give the corresponding calculations of the endurance strength 
distribution parameters for purposes of completeness and comparison. 
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Table 10 summarizes the endurance strength results used to prepare 
Fig. 13. 

2.2 Phase II Research 

The Phase II experimental research objective was to obtain 
distributional S-N and Goodman diagram data with specimens of the 
same steel as for Phase I but having a theoretical stress concentration 
factor of 2.34 instead of 1.42 for Phase I. 

2.2.1 Geometry and Hardness of Research Specimens 

Fig. 7 shows the geometry of the new specimens. The 
manufacturing processes were carefully controlled, to assure metallur- 
gical and strength similarity to Phase I specimens. Accordingly 
thirty-five (35) 9" long ungrooved specimens and thirty-five (35) 9" 
long grooved specimens for use in tensile tests, and one-thousand- 
fifty (1,050) 6" long specimens for use in fatigue testing were obtained. 

Measurements were made of the ungrooved and grooved 
specimens to verify the contractor's ability to meet machining tolerance 
and hardness requirements. The first set of tensile test specimens 
were not acceptable when the hardness and base diameters were found 
not to be within specifications, and the yield and ultimate strengths 
of the ungrooved specimens were significantly lower than those for the 
Phase I specimens . 

To assure that Phase I and Phase II specimens, made of 
two different lots of AISI 4340 steel, were similar except for the 
change in groove size and its effect on strength, a second set of 
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tensile test specimens was obtained. The specimen diameters and radii 
were measured with an Optical Comparator, at 20 magnification, located 
at the Arizona Gear Co. in Tucson. The results are given in Tables 11 
and 12. The surface finish at the base of the groove of the grooved 
specimens was observed through a microscope and visually compared with 
Johansen's precision gage blocks made by Pratt and Whitney because of 
the inaccessibility of the base of the groove to a Profilometer. 
Surface hardness measurements were made in the University's Metallur- 
gival Laboratory with a Wilson Rockwell C Hardness Tester using a 
150-K^ load and the Braile indenter. The results are given in Tables 
11 and 12. In addition, hardness tests were made on interior cross 
sectional areas of three grooved and three ungrooved specimens to 
determine the uniformity of hardness throughout each specimen. 

Standard precautions were taken to assure that the original hardness 
was not altered during sectioning. The location of test sites for the 
interior hardness measurements are shown in Figs. 17 and 18. The 
results of the interior hardness measurements are given on Tables 13 
and 14. 

Applying a 3o analysis to the surface hardness data in 
Tables 11 and 12 indicates that 99.73% of the population would lie 
between the limits of 36 R £ and 39 R £ for the grooved specimens, and 
between 35 R^ and 39 R c for the ungrooved specimens, hence within the 
specifications of R^ 35/40. The hardness readings in the sections 
given in Tables 13 and 14, show a preponderance of values between 
35.5 R c and 37.5 R £ with a high degree of uniformity. Since the 
hardness at the surface and in the interior of the specimens met the 
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specification requirements of 35/40, it was concluded that the 
specimens were given a proper heat treatment. 

2.2.2 Strength Characteristics of Research Specimens 

The ungrooved and grooved specimens were subjected to 
tensile loads to determine their yield, ultimate and breaking strengths. 
The tensile pull tests were performed at the Hughes Aircraft Company, 
Tucson, Arizona on a 60,000 lb. Tinius Olsen test machine. The machine 
was calibrated a short time before the tests and was considered to be 
in a fully operational condition. 

The thirty-five (35) ungrooved specimens were tested 
with an extensometer attached to each specimen which provided elong- 
ation input to a load-elongation pen recorder. The load was con- 
currently displayed on a 30-in. diameter dial segmented into to psi 
intervals on the 0 - 60,000 psi scale. The elongation of a two-inch 
gage length was measured with a micrometer after the specimen failed, 
and the percent elongation was calculated. 

As the specimen was placed under tension at a constant 
rate of elongation, observers watched the dial in an effort to identify 
the yield, ultimate and breaking loads. Upon reaching the ultimate 
point, the specimen was unloaded, the extensometer was removed, and 
the specimen diameter was measured with a micrometer to determine any 
reduction in cross-sectional area. The load was reapplied and the 
specimen was stressed to the breaking point. The yield load was 
identified on the recording at the 2% elongation point and the ultimate 
load was identified by the maximum load recorded on the chart. The 
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load at the breaking point, determined by visually observing the moving 
pointer, was manually recorded. The yield, ultimate, and breaking loads 
thus obtained are given in Table 15 and the percent elongation data 
and results are given in Table 16. 

After completion of the tensile test, the diameter of 
the ungrooved specimens were remeasured on the Optical Comparator to 
determine the area to be used in calculating the breaking strength. 

The final diameters and the elongation measurements are given in Tables 
15 and 16. The yield and ultimate strengths were calculated using 
the original (pre-test) specimen diameters; whereas, the breaking 
strength, was calculated on the basis of the final (reduced by elongation) 
diameters. The results are given on Table 15 and are summarized at the 
end of Table 15. 

During the testing of the grooved specimens the 
extensometer could not be used and the dial readings were visually 
observed and manually recorded. During the test of the first specimen 
it was observed that (1) the yield load could not be positively ident- 
ified, (2) there was little time lag and reduction in load between the 
ultimate and the breaking loads, and (3) the actual fracture of the 
specimen was accompanied by a loud shock wave throughout the testing 
laboratory. A decision was made thereafter to unload the specimens 
after the ultimate load was reached. Thus only the ultimate loads could 
be obtained, and are given in Table 17. Any change in specimen diameter 
at the ultimate load could not be measured accurately enough, hence 
it was decided to use the original (pre-test) area to calculate the 
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ultimate strength. The results are given in Table 18, and a 

\ 

summary thereof is given at the end of Table 17. 

2.2.3 Analysis of Tensile Test Results 

Before making the final decision regarding the accept- 
ability of the tensile test specimens as a basis for procuring fatigue 
test specimens, a review was made of the strength parameters obtained 
from test of the Phase I and Phase II specimens. The strength para- 
meters of the Phase I tensile test specimens, extracted from our 
previous report [3, Tables 7 and 8], and the results of the Phase II 
tests just discussed are listed in Table 18. It is noted that the 
mean yield and ultimate strengths obtained from Phase II tensile tests 
of ungrooved specimens are lower than the yield and ultimate strengths 
of the Phase I specimens. The standard deviations are different but 
within 4 percent of the respective means. The ultimate strength for 
the grooved specimens of Phase II was higher than for Phase I. This 
should be expected because Phase II specimens have a smaller groove 
radius which results in a higher static ultimate strength; consequently, 
this does not provide a basis for believing that the Phase II and 
Phase I specimens are different. 

The student "t" test was used to perform a comparison 
of the sample means of the ultimate strengths of the Phase I and 
Phase II tensile test specimens [7, pp. 193-194]. This method was 
selected because it is widely accepted to be valid for small sample 
sizes, and the Phase I data means were based on a sample size of 10. 

The results of the statistical tests and the corresponding critical 



values at the 0.05 level of significance are as follows: 



t - Statistic 

t - Statistic 



Critical Value 

Ungrooved Specimens 

19.65 

2.02 

Grooved Specimens 

13.10 

2.02 


The t-statistic is larger than the t-critical value which indicates 
that the difference between the means of the ultimate strength is 
statistically significant. The Phase I ungrooved specimens are 
apparently the stronger. 

The sample variances of the ultimate strengths were 
compared using the F-test [8, pp. 167-172]. The results, at the 0.05 
level of significance, are as follows: 

F - Statistic F - Statistic 

Critical Value 

Ungrooved Specimens 3.00 2.84 

Grooved Specimens 1.08 2.85 

The F-test for the variance further indicated a significant difference 
between the Phase I and Phase II ungrooved specimens but not between 
the grooved specimens. 

The results of the statistical analysis do not provide 
a logical basis for accepting or rejecting the Phase II specimens. 

Thus it became necessary to use different criteria to determine the 
acceptability of the tensile test specimens as a basis for ordering 
fatigue test specimens. 
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A metallurgical analysis of a specimen was obtained to 
determine its composition as AISI 4340 steel and its compliance with 
MIL-S-5000B. The analysis performed by Magnaflux Corporation, Materials 
Testing Laboratories confirmed that the sample met all requirements for 
the chemical composition. 

A review of the data from the physical measurements 
given in Tables 11 thru 14 confirmed that the test specimens met all 
specification requirements. The hardness values, which are considered 
to be most critical are well within specification requirements for 
surface hardness, and the additional hardness measurements on interior 
cross sectional surfaces also confirmed the uniformity of the hardness. 
All factors considered, the decision was made to proceed with the pro- 
curement of fatigue test specimens for Phase II. 

2.2.4 Recalibration of NASA Complex-Fatigue Research Machines 
2. 2. 4.1 Requirements 

Machine No. 2 required recalibration since it 
had been modified by the installation of spherical bearings at 
the Flex Couplings. Machines 1 and 3 were also recalibrated 
so that the strain gages of all machines would be verified to 
be functioning properly, and to revise the calibration co- 
efficients if the calibration results so indicated. The speed 
of each machine was also measured to assure its constancy and 
to determine if any change occurred due to wear in each machine 


and its drive motor. 
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2. 2. 4. 2 Procedures and Results 

The calibration method and procedures used 
were the same as those described in NASA CR - 72839 [4, pp. 
30-33] . The bending calibration was accomplished in two 
phases. First an out -of -machine bending calibration was done 
using the setup shown in Fig. 19. This phase verified that 
the bending bridge strain gages, shown in Fig. 3, mounted on 
the toolholder arm of each machine and the strain gage mounted 
in the groove of the calibration specimen were functioning 
properly and provided a relationship between the actual and 
the apparent strain in the toolholder strain gages. This was 
based on the facts that the actual /apparent stress ratio was 
close to the one it should be, and the plot of strain from the 
gage bridge versus applied weight was linear. 

The test setup shown in Fig. 20 was used to 
determine the calibration coefficients for the torque bridge, 
shown in Fig. 4, and for any interaction between torque and 
bending bridges . 

Next an in-machine, quasi-dynamic calibration 
was performed which provided the calibration coefficient needed 
to calculate nominal bending stress in the specimen groove 
from the apparent stress at the toolholder. 

Each machine was carefully calibrated with 
observations repeated a minimum of six times at each test 
point. The relation between the calibration variables, in 
each case, proved to be linear with the functional relationship 
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beginning at the origin. Thus the slope of each regression 
line for the calibration data completely defined the function 
and provided the calibration coefficients listed in Table 19. 

The rotating speed of each machine was 
determined using a tachometer strobe light which contains an 
internal oscillator and 60-cycle calibration. The results 
confirmed proper operation of each induction motor in providing 
a constant speed drive to each research machine. It was also 
found that the speed of each machine was independent of the 
bending and torque loads. The speed of each machine is listed in 
Table 19. The last calibration coefficients of the previous 
calibration was designated as Mode 4. The coefficients 
designated as Mode 5 apply to all data from June 1, 1971 until 
the next calibration. 

2.2.5 Ungrooved Specimens Fatigue Research 

Thirty-five Phase I specimens were machined down as 
shown in Fig. 21 so as not to provide any stress concentration at all. 
These specimens were tested in the Ann Arbor (R. R. Moore type) 
rotating beam fatigue research machine in our Reliability Research 
Laboratory to determine the endurance strength of ungrooved specimens. 

An optical comparator was used to determine the point 
of minimum diameter of each specimen and to measure the diameter. It 
was found that the specimens had a mean diameter of 0.3151 in., a 
standard deviation of 0.0004 in., and a range of 0.3143 in. to 0.3159 


in. 



A computer program was written and executed on the 
PDP-8 Computer to calculate the bending moments and pan loads required 
to test ungrooved specimens at target stress levels and desired 
specimens diameters.* The program was based on the following cali- 
bration equation for the Ann Arbor machine: 


M = 0.267 + 4.09L 


where 


(3) 


M = bending moment at the test section 
L = total effective load = pan weight (P) + 8.625 lb. 
Therefore the pan weight, P, to obtain a bending moment, M, is given 
by 


P 


(M + 0.267) 
4.09 


8.625 


(4) 


where 


M 


tt D x Stress 
32 


(5) 


and D is the test section diameter of the specimen to be tested next. 

In view of the range of specimen diameters, the 
computer program was run to calculate pan loads at diameter intervals 
of 0.0005 in. from 0.3140 in. to 0.3160 in. for each stress level 
desired. With interpolation, the required pan weights could be 
determined to the nearest +_ 0.10 lb. 

The staircase method of determining the endurance 
strength was used with a stress increment of 1,900 psi and a life of 
2.5 x 10^ cycles. A specimen was selected at random, its diameter 
obtained from the table of specimen diameteis, and the appropriate pan 


* This program is given in Appendix F. 
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weight as determined from the PDP program printout was used to apply 
the desired stress level. 

Thirty-eight valid runs were completed, as shown in 
the staircase plot of Fig. 22, with 15 successes and 13 failures. 

These data were then used to calculate the mean and standard deviation 
of the endurance strength distribution at 2.5 x 10^ cycles of life for 
the AISI 4340 steel, R c 35/40, ungrooved specimens subjected to an 
alternating bending stress and a constant shear stress with a shear 
ratio of r g = ~, as shown in Table 20. The mean endurance strength 

was found to be 80,725 psi and the standard deviation 3,040 psi. In 
comparison, the published endurance strength is estimated to be approx- 
imately 89,000. psi for polished specimens and 81,000 psi for machined 
specimens [9, pp. 160-172], The mean endurance strength as determined 
by this test program is very close to the published endurance strength 
of polished specimens. Thus the results of the tests are reasonable 
and compatible with current fatigue failure theory. 

When the decision was made to procure a second group 
of tensile test specimens, it was decided to obtain another group of 
ungrooved specimens and repeat the above tests to provide another 
basis for accepting the Phase II specimens as having physical prop- 
erties close to those of Phase I specimens. A number of Ann Arbor 
research machine outages had been encountered because the bending 
moments at higher stress levels approached the limit of the machine. 

Thus, it was decided to reduce the diameters of the new specimens to 
a nominal value of 0.2500 in. Thiry-five Phase II research specimens 
were obtained as per Fig. 21. The diameter of the specimens were 
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measured at the minimum point and the mean diameter was found to be 
0.2504 in. with a standard deviation of 0.0004 in. A specimen was 
randomly selected and the appropriate pan weight was determined from 
the PDP Computer program as before. The endurance life was again 
taken to be 2.5 x 10^ cycles and the staircase stress increment 1,900 
psi . 

Thirty -four useful data points consisting of 18 
successes and 16 failures were obtained, as shown in the staircase 
plot of Fig. 23. The 16 failures were used in the calculations given 
in Table 21. The endurance strength distribution parameters of the 
0.2500 diameter, ungrooved. Phase II research specimens subjected to 
an alternating bending stress and no shear stress, with a stress ratio 
r g = “, were found to be: Mean = 80,230 psi, and standard deviation 

= 1,425 psi. As the mean endurance strength was identical for both 
Phase I and Phase II steel specimens, the decision to continue using 
this steel and the manufacturer of the research specimens for the 
Phase II research was upheld. 

2.2.6 Grooved Specimens Fatigue Research 
2. 2. 6.1 Cycles-to-Failure Tests 

The Phase II test program using the new 
specimens, grooved to provide a theoretical stress concentration 
factor of 2.34, was initiated at the stress ratio of infinity. 
Using the von Mises-Hencky failure theory, the stress ratio , 

r s , for the loading provided by the research machines used in 
this research, is defined as 
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where 



( 8 ) 


S = alternating reversed bending stress due to the bending 
moment . 

t = mean shear stress due to torque, 
xym n 

Cycles-to-failure tests at the stress ratio of infinity 
were conducted at mean nominal alternating stress levels of 108,900, 
92,100, 73,600, 49,400, and 39,300 psi. Five stress levels are chosen 
to obtain data over the finite life range for the preparation of the 
corresponding S-N diagram. 

Upon completion of testing at the stress ratio of 
infinity, a decision was made to conduct all fatigue tests at a 
specific ratio on one test machine. Therefore, tests were initiated 
and completed at stress ratios of 1.06, 0.40, 0.25, and 0.15. Fewer 
stress levels were run as shown in Table 19 for stress ratios lower 
than » because of yield stress limitations in shear. 

The sample size to obtain the cycles-to-failure 
distribution for each alternating stress level and stress ratio com- 
bination was increased to 35, as recommended in the previous report 
[5, p. 95]. The actual alternating bending stress, the shear stress, 
the normal mean stress, and the stress ratio for each specimen in the 
sample were calculated by computer program STRESS. In addition, the 
mean and standard deviation of the achieved stresses for each sample 
were calculated, and are given in Table 22. 
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Program STRESS was updated for Phase II to incorporate 
current calibration constants for each machine to reduce the test data, 
and add rpm values for each machine. The computer printout includes 
a listing of the cycles to failure for each specimen in the sample. 

The updated program in Fortran language is given in Appendix B. 

Individual cycles-to-failure data were used as input 
data for program CYTOFR, as discussed in the previous report to NASA 
[5], to calculate the cycles-to-failure distribution parameters for 
the normal and log normal distributions and perform goodness-of-fit 
tests. Program CYTOFR was further updated to incorporate Cal-Comp 
graph and plot subroutines. Subroutine GRAPH constructs a histogram 
of the cycles to failure based on the failures per cell determined by 
the Chi-Squared test, and superimposes the distribution curve from the 
parameters computed by the main CYTOFR program on each histogram. The 
Chi-Squared test modification incorporated the automatic combining of 
cells at the tails of the distribution when the end cells do not 
contain at least five failure data points. The updated program in 
extended Fortran language is given in Appendix C. 

Program CYTOFR was run for each alternating bending 
stress level at which specimens were tested at the stress ratios of 
infinity, 1.06, 0.40, 0.25 and 0.15. Typical examples of the histograms 
and curves for the normal and log^ normal distributions are given in 
Figs. 24 and 25. The computed distribution parameters and values of 
the goodness-of-fit tests are summarized in Table 23 for the normal 
distribution, and in Table 24 for the log g normal distribution. 
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The goodness-of-fit test parameters listed in these 
tables were reviewed to determine if the applicable distribution of 
the cycles-to-failure data for Phase II specimens differed significantly 
from that for Phase I specimens. The following observations were 
made: 

1 . The values of the coefficients of skewness and kurtosis 
for the normal and the lognormal distributions were 
approximately the same and provided no preference for 
either distribution. 

2. The K-S goodness-of-fit test does not reject either 
distribution at the 0.05 level of significance. The 
largest D value was 0.201 for the normal distribution 
and 0.160 for the lognormal distribution; both were less 
then the critical D value of 0.224 at the 0.05 level 

of significance for a sample size of 35. Furthermore, 
there was essentially no difference between the normal 
and lognormal D values at any stress level. 

3. The Chi -Squared test proved to be more discriminating 
than the K-S test. The Chi-Squared value and the 
appropriate degrees of freedom (d.o.f.) for each stress 
level were compared with the critical value of 3.841 
for 1 d.o.f. and of 5.991 for 2 d.o.f. . The normal 
distribution was rejected in three out of seventeen 
samples; and the lognormal distribution was rejected 

in only two out of seventeen samples. One test each 
for normality and lognormality was inapplicable because 
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the data points were contained in only three cells 
resulting in zero degrees of freedom. 

4. It was concluded from these observations that there is 
a preference for the lognormal distribution over the 
normal distribution when working with cycles-to-failure 
data. 

The log g normal distribution parameters from Table 24 
were used to construct the S-N diagrams for the Phase II 
specimens. These S-N diagrams are given in Figs. 26 thru 
30. 


2. 2. 6. 2 Endurance Tests 

Endurance runs for Phase II specimens were 
completed using the "staircase" method at alternating bending 
to mean shear stress ratios of », 1.06, 0.40, 0.25 and 0.15. 
The staircase plots are given in Figs. 31 thru 35. The 
endurance strength distribution parameters were calculated in 
Tables 25 thru 29, and are summarized in Table 30. These 
parameters were used to complete the S-N diagrams of Figs. 26 
thru 30, and to construct the distributional Goodman diagram 
of Fig. 36. 
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3. THEORETICAL RESEARCH 

3.1 Generation of Finite Life Distributional Goodman Diagrams for 
Reliability Prediction 

A methodology for developing finite life distributional Goodman 
strength surfaces from cycles-to-failure distributions at specified 
alternating stress levels has been developed by Kececioglu and 
Guerrieri [10]. The Goodman strength surface shows the combinations 
of alternating bending stress and mean shear stress allowable to the 
design engineer. This study also investigated the applicability of the 
distortion energy and the maximum shear stress failure theories to 
determine which provided better correlation with the experimental 
data generated during Phase I. The finite life Goodman surface, 
developed using from two to five calculated strength distributions at 
specified stress ratios, can be used to construct distributions at any 
desired stress ratio and applied to probabilistic design. 

It is also found that the von Mises-Hencky ellipse effectively 
models Goodman diagrams for life greater than 10,000 cycles. However, 
when the equation for the von Mises-Hencky ellipse was modified from 

( S a /S e ) 2 * (V S u ) 2 * 1 < 9) 

to 

(V S e) 3 * (V S uV ■ 1 < 10 > 

and the values of a were determined from plots of finite life Goodman 
diagrams, the values ranged from 1.91 at 200,000 cycyles to 2.38 at 



40,000 cycles. This is an area that requires further study in search 
of a general equation, or applicable values of a, valid over all 
ranges of cycles to failure. 

3.2 The Weibull Distribution as A Description of Fatigue Life 
3.2.1. Introduction 

Since the Weibull distribution was introduced in 1949, 
it has gained wide acceptance as an extreme value distribution [11] . 

It has been used extensively in such areas as bearing fatigue data 
analysis [12], and the prediction of the failure of automotive parts. 

The Cal -Comp plots of normal and lognormal distributions 
did not fit the cycles-to-failure data very well for some sets of 
data. It was suggested that the three-parameter Weibull distribution 
may provide a better fit. 

The theory of the Weibull distribution and its manual 
application to experimental data are described in the existing 
literature [11], [12]. This study developed and validated a computer 
program which computes the parameters of the three-parameter Weibull 
distribution for a set of cycles-to-failure data, performs goodness- 
of-fit tests, and calculates the cycles to failure for 0.90 and 
0.99 reliability with 90% confidence. The program was tested using 
the data generated under Phase II and the results were analyzed. 
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3.2.2 Development 

The computer program development was patterned after the 
method described by Lochner [12], which used Weibull probability paper 
and manual calculations. The foundation is the generalized Weibull 
frequency distribution shown in Hahn and Shapiro [11, p. 110] 

f(N; 8, n, Y) = exp [ - (~^) 6 ] * (11) 

N >Y* -~<y < o, 8 > o, n>o, 

where 

N = cycles to failure 

8 = shape parameter or Weibull slope 

n = scale parameter 

y = location parameter . 

From the basic definition of reliability 


R 


{OO 

f (N) dN, 


the relationship between reliability and fatigue life N is 


R 



( 12 ) 


where y, n, and 8 are constants to be determined by the analysis of 
test data. The fraction failed, or unreliability Q, for N cycles 


is given by 




Equation (13) gives the probability of failure in N cycles or less, and 
is the cumulative distribution function, F(N). It can be transformed 
into a linear form by taking natural logarithms as follows: 



and 

In In [ j - _ - 1 F(N) ] = 3 In (N - y) - 6 In n. 

Letting 

y ■ ln ln h • 

and 


x = In (N - y), 


Eq. 14 becomes 


(14) 

(15) 

(16) 


y = 5 x + constant, 

or a transformed linear function. 

Weibull probability paper has been prepared with log log versus 

log scales so that the plot of y versus x of data would be a straight 

-(—) 0 

line with a slope of 6. When (N-y) = n, F(t) = 1 - e \n J =1 

- e * = 0.632, thus the value of (N-y) at which F(N) = 0.632 is an 
estimate of n- The location, parameter, y, is the minimum life point 
that provides the best approximation to linearity between x and y. 



3.2.3 Weibull Computer Program 

The basic FORTRAN computer program to determine the 
estimates of the Weibull parameters for cycles of life for specified 
levels of reliability was provided by Mr. Thomas C. Stansberry, Delco 
Radio Division, General Motors Corporation. His program was adapted 
to The University of Arizona CDC 6400 Computer and was updated to 
include subroutines for the Chi-Squared and Kolmogorov-Smimov goodness 
of-fit tests. Program WEIBULL is given in Appendix d. 

The first data card contains the sample size and the 
minimum life increment for use in linearizing the x-y relationship. 
Subsequent data cards (one for each specimen) contain the cycles -to- 
failure information. The first operation performed by the computer is 
to establish an ordered array of the cycles to failure and the 
corresponding median ranks. The computer calculates the median ranks, 
y i = In In q p (N ' y) > and x.^ = In (1L - y fc ) . Where i = 1, 2, ---, 

n and y^ = minimum life increment (1, 2, t k) such that y^ < . As 

the array of y^ and x^ is computed for different y^ the method of 
least squares is used to determine the degree of linearity. This 
operation is iterated with y^ being increased in increments until the 
best fit straight line is obtained. At that time the computer 
records the estimates of y, 3, and n. It then calculates the one 
percent failing life, the ten percent failing life, and the 50 percent 
failing life with the associated 90 percent confidence limits. 

Upon completion of the calculations, the program calls 
the K-S and Chi-Squared test subroutines, in turn, to provide a 
measure of goodness of the fit of the estimated Weibull distribution 
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with the data. The K-S subroutine, "DTEST", applies the Kolmogorov- 
Smirnov goodness-of-fit test and prints differences, D, for each 
, failure time in ascending order. Analysis of the K-S test is done by 
comparing the largest in absolute D value, with its critical value 
obtained from a D value table. 

The Chi-Squared test requires the subdivision of the cycles-to- 
failure data into a number of cells, k, determined by Sturges' rule 
[13] 

k =1+3.3 log 10 (n) , (17) 

where n is the sample size of the data. However, 
analysis of the results requires at least five data points in each 
cell. The use of Sturges' rule for data with a sample size of 35, 
results in six cells of equal width. Consequently, when the data is 
grouped into these six cells the cells at each end usually end up 
with fewer than five data points. If the two end cells are combined 
to provide five or more data points, the number of filled cells reduce 
to as few as four. Since the distribution being tested is the three 
parameter (r=3) Weibull, the degrees of freedom (k-r-1) requires that 
the number of cells, k, be at least five in order to have at least one 
degree of freedom. This Chi-Squared test was applied to samples of 
Phase II cycles-to-failure data. It was found that six of the 
twelve tests resulted in only four filled cells. Thus, there was zero 
degrees of freedom and the Chi-Squared test was not useable. Thus, 
it appears that the sample size will have to be increased still further 
if the standard Chi-Squared goodness-of-fit test is to be used. 
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To circumvent this problem variable cell widths were 
used. The technique described by Hahn and Shapiro [11, pp. 302-308] 
called for the calculation of cell widths to provide equal number of 
observations in each cell. A modification to the subroutine was made 
dividing the range of the data so that each cell contains exactly five 
cycles-to-failure data. For our sample size of 35, this provides 
seven cells. This subroutine was run for the same 12 sets of data. 
The expected frequency for the seventh cell was always less than five, 
which according to accepted practices invalidate the test. It was 
observed, however, that as long as the expected frequency was equal 
to or greater than two, the Chi -Squared value could be calculated. 

This observation was confirmed by a Monte Carlo Simulation of 1,000 
runs from which it was concluded that Chi -Squared errors resulting 
from expected frequencies between two and five are insignificant. 
Nevertheless, a final modification was made to the subroutine using 
variable cell widths, but combining adjoining end cells to insure that 
the expected number of observations per cell is equal to or greater 
than five. When the same cycles-to-failure data was rerun to apply 
this Chi-Squared test subroutine, all 12 tests resulted in six useable 
cells thus providing two degrees of freedom. 

3.2.4 Results 

The operation and accuracy of the computer program were 
verified by using the same input data used by Lochner [12]. Identical 
estimates were obtained for each parameter to the degree of accuracy 
obtainable from probability paper plots. The computer program 



provided parameter estimates to five place accuracy and used this 
accuracy in subsequent calculations. In Table 31 the Weibull distri- 
bution parameters, and the K-S and Chi -Squared goodness -of-f it results 
are given. A sample Cal-Comp plot of the Weibull distribution is given 
in Fig. 37. 


The accuracy of the D-test subroutine for the Kolmogorov- 
Smirnov goodness-of-fit test was confirmed by a desk calculator. The 
maximum D values found by applying the subroutines to the cycles-to- 
failure data are listed in Table 31. These results show that, at the 

0.05 level of significance, the K-S test does not reject the Weibull 
distribution in all of the 17 tests. Thus, the Weibull distribution 
can safely be used to approximate distributions of cycles-to-failure 
data at specified stress levels. 

The Chi-Squared test values determined by the WEIBULL 
subroutines with variable cell widths are also given in Table 31. 

Note that the variable cell width analysis rejects 5 out of 17 tests. 

The conclusions drawn from the above analysis are listed 

below: 

1. Based on the K-S test results of not rejecting any of the 
samples, the three-parameter Weibull distribution may 
describe fatigue cycles-to-failure data. 

2. Based on the Chi -Squared test results of 5 rejections out 
of 17 samples, the Weibull distribution may not be 
considered generally acceptable for cycles-to-failure data 
consisting of 35 data points. Further in 14 out 16 cases 
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the Chi -Squared value for the Weibull is greater than for 
the lognormal. 

3. Based on the previous two conclusions and the results in 
Tables 24 and 31, the lognormal distribution appears to 
represent the cycles-to-failure data of the Phase II 
research best. 

3.3 Reliability of Components Subjected to Cumulative Fatigue 
3.3.1 Introduction 

In cumulative fatigue of most concern to design 
engineers is the mathematical relationship between the number of load 
cycles at various alternating stress levels applied to a component, 
the S-N diagram results, and the survival life of the component under 
these conditions. After such a relationship is determined a method 
needs to be developed to predict the reliability of a component sub- 
jected to a specified history of cumulative fatigue stresses. The 
objective of this study was to review the published cumulative fatigue 
theories, and to discover or come up with methods for making reliability 
predictions . 
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3.3.2 Literature Search 

The literature search revealed attempts to describe the 
degree of cumulative damage in expressions involving transfer of 
energy or mass, with damage often described by a crack parameter and 
interpreted by Osgood [14] as a change in the state of energy in the 
immediately adjacent volumes of material. The primary difficulty 
with these methods is their complexity and highly approximate nature. 

The simplest and most widely used cumulative damage rule 
is Miner's rule [15], based on the assumption that cumulative damage 
under cyclic stressing is related to the net work absorbed by the 
specimen. That is, the number of stress cycles applied, expressed as 
a percentage of the number of cycles of life at the given alter- 
nating stress level, is the proportion of useful life expended. 
Therefore, the specimen should fail when the total damage reaches 1.00, 
or 

m n. 

* IT” = 1 < 18 > 

1=1 i 

where 

n^, n^, . . . , = cycles of operation at each applied 

alternating stress level. 

N, , N_ . . . , N = cycles of life at each stress level. Miner's 

12 m 
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rule is accepted as providing a good, conservative first approximation 
for engineers in preliminary design, but fails to account for the 
effects of overstressing or understressing in the early cycles or for 
loading sequence. 

An expression was developed by Corten and 

Dolan [16] to model the hypothesis that fatigue damage in terms of the 
nucleation of submicroscopic voids which develop into cracks, is a 
function of damage nuclei and the rate of damage propagation. Damage, 
which was represented as a power function of the number of cycles, 
was summed for a loading sequence consisting of repeated blocks of 
cycles alternating between two stress amplitudes. The functional 
relationship developed is 


N 

g 


°1 + a 2 ( S 


N 1 

S 2 .d f S 3 .d 

•) ♦ °3 (“§-) + 


1 


'1 


, n >d 

. . + a (•= — •) 
n o 2 


where 


(19) 


N = total number of cycles of stress to failure for an 
g 

incremental stress amplitude history, 

= number of cycles at the highest stress level, S^, before 
failure , 

a^, a^, • • . , a = ratio of the number of cycles applied at 

stress levels S^, S 2 , .... to the 
total cycles applied, 

Sj > > ... > = various alternating stress levels or 

amplitudes applied , 

d = inverse slope of the linear portion of the S-N diagram. 
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The Corten-Dolan method appears to give a better 
correlation with existing data than Miner's rule; however, it still 
has the deficiencies that the value of d cannot be determined with a 
reasonable accuracy, and the equation is based on a deterministic 
rather than a distributional S-N diagram. 

The NERVA program [17] approaches the cumulative fatigue 
problem by revising Miner's rule as follows: 

m n. 



where 

y = a normally distributed variable with a mean, y, and a 

standard deviation, o . 

Y 

The experimental values of y have been found to range 
between 0.18 and 23.0, depending upon the material, the test conditions 

and the order of loading history. It was observed that a low to high 

\ 

loading sequence (s^ < < s^..) resulted in high y values (y > 1). 

A high to low loading sequence (s^ > s^ > s^ . ..) gave low y values 
(y < 1). For a loading history with representative high, low, and 
medium stress levels in random order the value of y appears to be 
close to unity. 

Sorensen [18] developed a general expression for the 
probability distribution of the damage rate in 

terms of the power spectrum of the random excitation. The method 
requires use of the single valued theoretical S-N diagram to determine 
distributional values. The analysis is logical and the results 


obtained from a numerical example appear to be reasonable. The 
method warrants further investigation using distributional S-N diagrams 
developed by this research. Serensen [19, pp. 33-43] studied fatigue 
damage accumulation under distributional service loading. A family 
of curves for the probability density function of the stress 
amplitude, the distribution of fatigue life for a given group of parts 
at a given type of loading, and Miner's rule are used to assess fatigue 
damage. The analysis and probabilistic calculations are logical but 
are subject to the limitations of Miner's rule. 

3.3.3 Proposed Methods to Calculating the Reliability of 

Components Subjected to Cumulative Fatigue at Sequenced 
Stress Levels. 

Use was made of the cumulative fatigue theories and the 
statistical nature of fatigue life to develop methodologies for the 
calculation of reliability. Two methods were developed. 

The first method makes use of the multiplication rule 
and conditional reliabilities or probabilities of survival. With 
this approach the first step is to calculate the probability, P of 
surviving the first cycles at stress level . The next step is 

to compute the probability, P , of surviving cycles at stress level 

given survival of cycles at stress level . Then the 

probability of survival for the sum of these cycles (N^ + N^) is the 
product of the individual probabilities (P^ • ^ 2 ^' 'This procedure 
would be continued for as many steps as necessary. 



The second method is called the method of equivalent 
reliabilities. Here the probability of surviving cycles at stress 
level S 1 is computed. Utilizing this information an equivalent cycle 
life, , at stress level S 2 is computed. The N 2 cycles at stress 
level S 2 are now added to NJ and the probability of surviving the Nj 
+ N 2 cycles at stress level S 2 is computed. This value is equivalent 
to the reliability associated with the survival of the cycles at 
stress level and the N 2 cycles at stress level S 2 * Again, this 
procedure could be continued as many times as needed to obtain the 
final reliability of a component subjected to cumulative fatigue at 
sequenced stress levels. 

The two methods will be evaluated next using life cycle 
data obtained from the experimental test program. Consider the 
following stress history applied to a steel shaft: 


Alternating Stress Level 

Cycles Run 

= 86,000 psi 

10,000 

S 2 = 96,000 psi 

1,000 

S 3 = 100,000 psi 

500 


The fatigue life of specimens tested at these mean 


stress levels were found to be as follows: 



40 


Alternating Fatigue Life 


Stress Level 

Mean Cycles 

Standard Deviation 


Los 10 

Cog jo Cycles 

= 86,000 psi 

4.715 

0.068 

S£ = 96,000 psi 

4.394 

0.052 

S_ =100,000 psi 

4.102 

0.073 


Consider first the method that uses the multiplication 
rule and conditional probabilities. With the stress levels increasing, 
the reliability for the first level is given by 


fco 


fj (Nj) dN 1 , 
' N 


( 21 ) 


where f^CN^) is the normal probability density function (pdf) for the 
log 10 cycles-to-failure at stress level . Transforming to the 
standard normal pdf variable 


, N-N , 

z = (— ), gives 

N 


R 


1 


(°° 

<f> (z) dz . 


1 z 


1 


In this example 


log 10 (10,000) -4.715 
0.068 


-10.5, 


( 22 ) 


which from standard normal distribution area tables yields a relia- 
bility, Rj, of essentially 1 from Eqs . (2l) and(22j. 



The conditional reliability for the second stress level 


is given by 


r 2 = J 


N W dN 2 
_2 


(23) 


f,(N 0 ) dN„ 


N 


1 


where ^^2^ * s t * ie c y c l es- to-failure at stress level S^. 

Transforming to the standard normal pdf variables results in 


where 


and 


R, 


<J> (z) dz 


(24) 


d> (z) dz 




log 1Q (11,000) -4.394 
0.052 


log (10.000) -4.394 
0.052 


= -6.75, 


= -7.54 


Insertion of these values of z^ and z^ into Eq. |24)yields a reliability 
R 2 , of 0.9 1q * 

Continuing in a similar manner for the third stress level 
the reliability is given by 


The symbol °- 9 io is used to represent the number of 0.9999999999. 
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or 


N, 


W dN 3 


R, * — 

3 '' 0o 


N, 


W dN s 


R. 


<!> (z) dz 


4> (z) dz 


where 


and 


log 10 (ll,500) -4.102 
0.073 


log 10 (ll,000) -4.102 
0.073 


The value for is found to be 0.9004. 


(25) 


(26) 


The reliability of the shaft for the mission length of 

11,500 cycles and the given alternating stress history is obtained by 

R = R. R 0 R,, which results in R = 0.90003. 

S 1 l o s 

The reliability was calculated with the stress 

history reversed; i.e. decreasing stress level. The resulting 

reliability was 0.9^. 
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The same problem is approached next utilizing the method 
of equivalent reliabilities. The increasing stress history is 
considered first. Calculating the z value for the first stress level 
gives 


log (10,000) -4.715 

z. = — = -10.5 . 

0.068 

Utilizing this value the equivalent number of cycles at 
stress level, S^, is found as follows: 

Z 1 = 2 2 = " 10 - 5 » 

where 

log 10 -4.394 

2 ? = » 

0.052 


so that 

log 10 = z' 2 (0.052) + 4.394. 

Solving for results in 7,050 cycles. This means that 
7,050 cycles run at stress level is equivalent to 10,000 cycles at 
stress level S^. Thus, to find the reliability associated with both 
stress histories, add the 1,000 cycles run at stress level to the 
equivalent number of cycles and calculate the new z value, or 

log 1() (8,050) -4.394 

Z 2 ~ 


0.052 


-9.38 . 
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In like manner the equivalent number of cycles at stress 
level is calculated as follows: 

z 2 = z' 5 = -9.38, 

_ 1Og 10 N 3 - 4 ' 102 

Z 3 ' 

° 0.073 


Log 1() =z^(0.073j+ 4.102, 

or 

= 2,610. 

Adding the 500 cycles run at stress level and 
computing the z value, the reliability is found to be 

N 3 = + 500 = 3,110, 



or 


R 


f(z) dz. 


z 


3 


R = 0.9 


15 ' 


(27) 


The method when applied to the stress history applied in 
decreasing order yields a z^ value of -8.36 and a corresponding 
reliability of 0.9 . 
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Comparison of the two methods indicates the following: 

1. The method utilizing conditional probabilities gives a substantially 
lower reliability with the increasing stress history than with the 
decreasing stress history. The large difference in these values is 
unexpected since it contradicts current theories of cumulative fatigue 
which ignore the effects of the order of stress history on the life 
of a component. 2. The method of equivalent reliabilities gives 
results which are consistant regardless of the order of application. 
There remains a need to experimentally verify the applicability of 
these methods. 

3.3.4 Conclusions 

Cumulative fatigue and its effect on component reliability 
are quite complicated and not completely understood areas. Several 
theories have been developed with most research devoted to understanding 
changes in material structure due to microscopic crack nucleation and 
propagation. Miner's rule [15] provides a rough approximation for 
determining the expected mean life of a component subjected to various 
stress levels. However, this rule does not provide a method for 
calculating component reliability. 

The equation developed by Corten and Dolan [16] provides 
an improvement of Miner's rule in that it correlates better with 
existing data. This method suffers from lack of accuracy in providing 
a deterministic value from distributional data. 

Work done under the NERVA program [17] to gain insight 

when stresses are applied in various sequences, has potential when 
combined with distributional fatigue life data. 
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The studies performed by Sorensen [18] and Serensen [19] 
approach the problem of fatigue using randomly applied or distributed 
loads resulting in distributed stress, which would constitute a 
generalized approach to cumulative fatigue. Their methods should be 
investigated further using our distributional fatigue data. 

The first method proposed here for calculating cumulative 
fatigue reliability combines the product rule and conditional prob- 
ability theory for sequentially applied levels of alternating stress. 

This method gives significantly different results for sequentially 
decreasing than increasing stress levels, namely O-^i anc ^ 0*9003 respect- 
ively. Thus this method does not appear to provide a general enough model. 

The second proposed method is that of equivalent relia- 
bilities where the equivalent number of cycles at each succeeding 
stress level is found for the reliability calculated at the preceding 
level. It includes consideration of fatigue damage at all stress 
levels and the effects of increasing versus decreasing stress history. 

The calculated reliabilities are essentially the same, thereby 
demonstrating consistency. This method is considered to provide a 
protentially valid general reliability model for cumulative fatigue. 



4. 


OVERALL CONCLUSIONS 


A methodology has been developed for designing specified 
reliabilities at optimum size and weight into rotating mechanical 
components subjected to fatigue under combined alternating 
bending stress and constant shear stress. 

Three research machines have been designed and fabricated which 
are capable of simultaneously applying desired levels of 
alternating bending moment and constant torque to rotating test 
specimens . 

Phase I of an experimental fatigue research program to verify the 
probabilistic design methodology was conducted with AISI 4340 
steel specimens grooved to provide a theoretical stress concen- 
tration factor of 1.42. The data obtained were reduced using 
three computer programs developed for the CDC 6400 computer. The 
results were used to construct a distributional Goodman strength 
diagram for 2.5 x 10^ cycles of life and distributional alter- 
nating bending stress versus cycles-to-failure (S-N) diagrams. 

The following specific conclusions were reached: 

3.1 The sample sizes of 12 and 18 were not sufficiently large 
for goodness-of-fit tests to determine whether the normal 
or the log-normal distribution provided a better fit to 
the data. The Chi -Squared test could not be used all the 
time and the Kolmogorov-Smirnov test did not reject either 
distribution. Phenomenological considerations, combined 
with goodness-of-fit test results, and coefficients of 
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skewness and kurtosis indicated that the log normal dis- 
tribution provided the better fit to cycles-to-failure 
data, while the normal distribution provided the better 
fit to endurance strength and finite life stress-to-failure 
data. 

3.2 Probabilistic S-N diagrams plotted for stress ratios of ®, 

3.5, 0.83 and 0.44 for Phase I data showed a linear re- 
lationship between the log of cycles-to-failure and the 
log of alternating bending stress levels for each stress 
ratio. 

3.3 There is significant reduction in mean fatigue life as 

the alternating bending strength is reduced from °° to 

0.44. For the alternating stress level of 70,000 psi 

the estimated cycles to failure, N^, are 200,000 for an 

r of <*>, N_ = 84,000 for r of 3.5, N_ = 80,000 for 

s f s f 

r of 0.83, and N_ = 51,000 for r of 0.44. 
s f * s 

3.4 There is a relatively high variability in the cycles-to- 
failure data with a coefficient of variation ranging from 
15% to 25%. However, there is sufficent consistency to 
provide the essential linearity in the probabilistic 
envelope of the S-N curve, as well as in the plot of the mean 
log cycles. 

3.5 The distributional Goodman strength diagram for Phase I 
results provides data directly useable by the designer for 
determinining the reliability of rotating components, for 
a service life of 2.5 x 10^ cycles, subjected to combined 
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alternating bending stress and constant shear stress at 
various stress ratios for the given material and geometry 
of the specimen tested. 

4. Phase II of the experimental fatigue research program was 
initiated and completed for specimens of identical specifications 
to the Phase I specimens, except for a different groove radius to 
provide a theoretical stress concentration factor of 2.34. 

Fatigue tests were accomplished for plotting S-N diagrams for the 
stress ratios of «>, 1.06, 0.40, 0.25 and 0.15. 

5. The computer programs were revised to incorporate the current 
calibration coefficients for the fatigue research machines, the 
application of a revised Chi-Squared goodness-of-fit test, and a 
subroutine for the plotting of a histogram and the estimated 
normal, lognormal, and Weibull distributions using a Cal-Comp 
plotter. 

6. The conclusions reached from the results of the computer programs 
are the following: 

6.1 Increasing the sample size to 35 allowed the use of the 
Chi -Squared goodness-of-fit test which proved to be more 
discriminating than the Kolmogorov-Smirnov test. The 
Chi-Squared test rejected the normal distribution three 
times out of sixteen samples at the 0.05 level of signifi- 
cance; whereas the lognormal distribution was rejected only 
two times out of sixteen. Thus, the previous tenuous 
preference for using the lognormal distribution for cycles- 
to-failure data was strengthened. For the Chi-Squared test 
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to be valid, it must have five or more observations in each 
cell; this requires that the end cells sometimes be combined. 
The existing methodology describes the tails of distributions 
with about 5% accuracy. A method to improve the accuracy 
of describing the tails of distributions for use in calcu- 
lating reliability is needed. 

6.2 Comparison of the S-N and Goodman diagrams for Phase I and 
Phase II research reveals that the endurance strength for 
2.5 x 10^ cycles of life and the finite fatigue life at a 
specified alternating stress level are significantly 
reduced when the stress concentration is increased by 
reducing the groove radius. The S-N diagrams show that 
the fatigue life for an alternating stress level of 80,000 
psi decreases from 75,000 cycles for Phase 1 specimens to 
9,000 cycles for Phase II specimens for a stress ratio of °°. 
In Phase I tests, an alternating stress level of 50,000 psi 
at the stress ratio of 0.44 approached the endurance level, 
and a mean cycles to failure could not be determined. For 
the Phase II tests at the stress ratio of 0.40, the 50,000 
psi alternating stress level approached the highest stress 
level that could be tested without suffering shear failures 
during machine set up, and at this stress level the mean 
cycles-to-failure was about 56,000 cycles. 

6.3 The cycles-to-failure results emphasize the effect of stress 
ratio on the fatigue life of specimens at specified levels 
of alternating stress and the maximum stress levels at 
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which fatigue tests can be conducted. The upper limits for 
the tests ranged from about 120,000 psi for alternating 
bending stress for the stress ratio of » to approximately 
32,000 psi for the stress ratio of 0.15. Examination of 

1 

the S-N diagrams show that there is a limited useable range 
for low stress ratios; consequently, the stress level and 
stress ratio become very critical in a component design 
using low stress ratios. 

7. Generation and analysis of finite life distributional Goodman 
diagrams reveal that the conventional deterministic Goodman diagram 
is extremently conservative. 

8. The FORTRAN Computer program to estimate the parameters of the 
three-parameter Weibull distribution which would represent cycles- 
to-failure data is providing good results. The program incorp- 
orates the K-S and Chi-Squared goodness-of-fit tests. The K-S 
test did not reject the Weibull in 17 tests, and the Chi- Squared 
test rejected the Weibull 5 times in 17 tests. Based on these 
test results, the Weibull can be considered a useable distribution 
for cycles-to-failure distributions, although the lognormal is 
favored by the Chi-Squared test. 

9. Cumulative fatigue is a complicated and not completely understood 
design area. Several theories have been developed with most 
research devoted to understanding changes in material structure 
due to microscopic crack nucleation and progression. Two methods 
were investigated in this research for calculating the reliability 
of a component sequentially subjected to a number of stress levels. 
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The method of equivalent reliabilities provides consistent 
component reliabilities regardless of the order of application 


of the different stress levels. 
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5. RECOMMENDATIONS 

1. The Phase II type experimental research program should be 
repeated for other geometries and materials, and a complete set 
of S-N and Goodman diagrams should be prepared. 

2. An analytical study should be made in search of a mathematical 
relationship between infinite life and finite life Goodman 
diagrams . 

3. An analytical study should be made in search of a mathematical 
model for the effect of varying stress concentrations and stress 
ratios in fatigue tests. 

4. Further study should be made of infinite and finite life Goodman 
strength diagrams to determine (1) a more accurate model than 
the ellipse and (2) if another failure theory, or combination of 
theories, is more valid for fatigue failure than the distortion 
energy theory. 

5. A study should be made with sample sizes larger than 35 to 
determine the effect of sample size on describing the tails of 
the statistical distributions. As part of the study an error 
analysis should be made comparing the predicted reliabilities 
using the lognormal and Weibull distributions for the cycles to 
failure data and for different sample sizes. 

6. A Phase III experimental research program should be conducted 

using ungrooved specimens so that the true effect of stress con- 
centration on fatigue life may be better quantified comparatively. 
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7. An experimental cumulative fatigue research program should be 

conducted by sequentially applying various levels of stresses for 
specified numbers of cycles and to specimens with different 
groove radii. The data thus obtained should be used to determine 
the validity of the equivalent reliabilities cumulative fatigue 
model . 
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Fig. 1 Schematic front view of combined reversed bending-steady torque 

fatigue-reliability research machine. 
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Fig. 2 Schematic top view of combined reversed bending- steady torque 

fatigue-reliability research machine. 
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Bending stress strain gage bridge. 
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Shear stress strain gage bridge. 
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Fig. 5 


Flow chart of steps in cycles -to- failure and stress-to-failure 
fatigue tests. 
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Fig 7 

r* , * 


Phase II research specimen geometry. 









Log 10 CYCLES TO FAILURE 

Cycles-to-failure distributions at the stress ratio of 
for AISI 4340 steel R c 35/40 Phase I grooved specimens 






Cycles-to-failure distributions at the stress ratio of 0.83 
for AISI 4340 steel R 35/40 Phase I grooved specimens. 
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Fig. 12 Endurance strength data obtained by the staircase method for stress ratio of 0.45 

for AISI 4340 steel, MIL-S-5000B, Condition C4, Rockwell C 35/40, with Phase I 
grooved specimens . 
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Fig. 14 Endurance strength data obtained by the staircase method for stress ratio of 

for AISI 4340 steel, MIL-S-5000B, Condition C4, Rockwell C 35/40, with Phase 
grooved specimens . 



O success If specimen survived 2.5 


74 



isd)| ‘uoi}38s asa} ui sssj^s Sutpuaq xbutuion 


Fig. 15 Endurance strength data obtained by the staircase method for stress ratio of 3.5 

for AISI 4340 steel, MIL-S-5000B, Condition C4, Rockwell C 35/40, with Phase I 
grooved specimens . 
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Fig. 16 Endurance strength data obtained by the staircase method for stress ratio of 1.0 

for AISI 4340 steel, MI L-S- 5000B, Condition C4, Rockwell C 35/40, with Phase I 
grooved specimens . 
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(a) The length AB used for hardness measurements. 


90° axis 


0°axis 


(b) Test pattern for section A 



Fig. 17 


Grooved specimen sections used for hardness measurements. 


Legend: 


C - center area 

I - intermediate area 

0 - outside area 1/16" from surface 


(a) Ungrooved specimen 



A (b) Grooved specimen 



Fig. 18 


Specimen sections and locations for hardness measurements. 



Torque clamp-screw device 



Bending stress bridge calibration setup. 
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Fig. 20 Torque (shear) stress bridge calibration setup. 




Fig. 21 Ungrooved endurance strength specimen. 




Fig. 22 Endurance strength data obtained by staircase method in Ann Arbor machine for the 

stress ratio of <*> for AISI 4340 steel R £ 35/40 Phase I research. Specimens with 
no groove and a test section diameter of 0.3150 in. 
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Fig. 23 Endurance strength data obtained by the staircase method in the Ann Arbor machine 

for the stress ratio of °° for AISI 4340 steel R £ 35/40 Phase II research. 
Specimens with no groove and a test section diameter of 0.2500 in. 








FREQUENCY/CLASS INTERVAL WIDTH 


NORMAL DISTRIBUTION PARAMETERS 



CYCLES-TO-FRILURE X 10 ’ 

MEAN VALUE: 3297.1 CYCLES 

STANOARO DEVIATION: 361-7 CYCLES 

KOLMOGOROV-SMIRNOV TEST: 0.094 
CHI-SQUARED TEST: 0.248 

SKEWNESS: -0.163 

KURTOSIS: 2-317 

FIG. 24 CYCLES-TO-FAILURE DISTRIBUTION OF 35 

GROOVED SPECIMENS FOR AN ALTERNATING 
STRESS LEVEL OF 108. 900 PS I AT 8 
STRESS RATIO OF INFINITY AND NOMINAL 
GROOVE DIAMETER OF 0-491 INCHES. 



FREQUENCY/CLRSS INTERVAL NIDTH 


LOG NORMAL DISTRIBUTION PARAMETERS 



CYCLES-TO-FAILURE 

MEAN VALUE: 8.095 CYCLES 

STANDARD DEVIATION: 0.112 CYCLES 

KOLMOGOROV-SMIRNOV TEST: 0.096 
CHI-SQUARED TEST: 2.844 

SKEWNESS: -0*374 

KURTOSIS: 2.408 

FIG. 25 CYCLES-TO-FAILURE DISTRIBUTION OF 35 

GROOVED SPECIMENS FOR AN ALTERNATING 
STRESS LEVEL 0F108.900 PSI AT A 
STRESS RATIO OF INFINITY AND NOMINAL 
GROOVE DIAMETER OF 0.491 INCHES. 
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Fig. 27 Cycles-to- failure distribution at the stress ratio of 1.06 

for AISI 4340 R 35/40 Phase II grooved specimens. 
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Cycles-to-failure distribution at the stress ratio of 0.40 
for AISI 4340 R c 35/40 Phase II grooved specimens. 




! scJ » *SS3yiS 9NllVNd3nV 3VNIW0N 


Fig. 29 Cycles-to-failure distribution at the stress ratio of 0.25 
for AISI 4340 R c 35/40 Phase II grooved specimens. 
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Number of specimen 

Pi„ 31 Endurance strength data obtained by the staircase method for Phase II grooved 

specimens of AISI 4340 steel Rq 35/ 40 for stress ratio of °°. 













Endurance strength data obtained by the staircase method for Phase II grooved 
specimens of AISI 4340 steel R 35/40 for stress ratio of 1.06. 










































Fig. 33 Endurance strength data obtained by the staircase method for Phase II grooved 

specimens of AISI 4340 steel R c 35/40 for stress ratio of 0.40. 













Endurance strength data obtained by the staircase method with Phase II grooved specimens 
of AISI 4340 steel R c 35/40 for stress ratio of 0.25. 













Endurance strength data obtained by the staircase method for Phase II grooved 
specimens of AISI 4340 steel R c 35/40 for stress ratio of 0.15. 
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FREQUENCY/CLASS INTERVAL WIDTH 


WE I BULL DISTRIBUTION PARAMETERS 



X 10 

KOLMGGOROV-SM I RNOV TEST: 0.217 

CHI-SQUflRED TEST: 4.236 

WE I BULL SLOPE ( BETA 3 : 3.853 

MINIMUM LIFE (GAMMA): 1999 

SCALE PARAMETER (ETA): 1306 


FIG. 37 CYCLES-TO-FAILURE DISTRIBUTION OF 35 GROOVED 

SPECIMENS FOR AN ALTERNATING STRESS LEVEL OF 
108,900 PSI AT A STRESS RATIO OF INFINITY" 

AND NOMINAL GROOVES DIAMETER OF 0.491 INCHES. 
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Table 1 Summary of nominal stresses applied to eighteen 

Phase I specimens for a nominal S = 65,000 psi 

and r e 0. 44. a 

s 



Mean 

Standard deviation 

Bending Stress (psi) , a = S 

xa a 

64,700 

1,300 

Shear Stress (psi) , t 

xym 

s 

84,900 

1,800 

_ . a 

Stress Ratio = 

T 

xym 

0.44 

0.01 


Table 2 Estimates for cycles- to-fai lure distribution 

parameters for S = 65,000 psi and r ~ 0.44. 


Distribution 

Normal 

Lognormal 

e 

Mean Cycles -to-Fai lure 

83,700 

11.278 

Standard Deviation 

26,350 

0.374 

Skewness 

0.315 

1.091 

Kurtosis 

2.408 

3.747 

K-S Test Maximum D Value 

0.066 

0.114 

Allowable D Value at the 
0.05 significance level 

! 

0.294 

0.294 
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Rounded off to nearest 100 psi. 

s m = t xym and <Ts m = /j ^Txym usin 8 the von Mises-Hencky failure theory. 


























Table 4 Normal distribution parameter estimates and the max D values for Phas 

cycles-to-failure results. 
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Table 6 Endurance strength calculations for the stress ratio of 

0. 45 using the staircase method at 2.5 x 10^ cycles 
for AISI 4340 steel R 35/40 Phase I grooved specimens. 


Nominal 

alternating 

s tress 



d = stress increment = 2,410 psi 

lowest stress level = 44,640 psi 
X = mean (estimate) 

X = X + d [A/N + 1/2] = 44,638 + 2,411 [ ii + I] 

_ ° 18 2 
X = 49,597 psi 2 49,600* psi 

s = standard deviation (estimate) 

s = 1.620 d[(NB-A 2 )/N 2 + 0.029] =. 1.260 (2,411) 

s = 3,681 psi s 3,700* psi 


[ 18(60) - (28) £ » 0 . 029 ] 
(IS ) 2 


* Sounded off to nearest 100 psi. 
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Table 7 . Endurance strength calculations for the stress ratio 

of « using the staircase method at 2.5 x 106 cycles 
for AISI 4340 steel R„ 35/40 Phase I grooved 
specimens. 



d = stress increment = 3,314 psi 
X q = lowest stress level =51,915 psi 
X = mean (estimate) 

X = X q + d[A/N + 1/2] =.51,915 + 3,314 [ l J + i ] 

X = 57,317 psi 2 57,300* psi 15 2 

s = standard deviation” (estimate) 

S = 1.620 d[(NB-A 2 )/N 2 + 0.029] * 1.620 (3,314) J 15(27) - (17) 2 + Q Q2Q 
s = 2,924 psi s 2,900* psi (15) 2 


* Rounded off to nearest 100 psi. 
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Endurance strength calculations for the stress ratio 
of 3.5 using the staircase method at 2.5 x 10^ cycles 
for AISI 4340 steel R c 35/40 Phase I grooved 
specimens . 


Alternating 

stress 

psi 



57,103 

3 

53,759 

2 

50,415 

1 

47,071 

0 



d = stress increment = 3,344 psi 

X = lowest stress level = 47,071 psi 

_o 

X = mean (estimate) 

X = X 4 d[A/N + 1/2] = 47,071 4 3,344 [M + I ] 
o J 18 2 

X 55,059 psi. — 55,100* psi 

s = standard deviation' (estimate) 

S = 1.620 d[(NB-A 2 )/N 2 4 0.029] = 1.620 (3,344) [ 18 ( 76 ) - (34) 2 + 0.029] 

s = 3,701 psi = 3,700* psi (18) ^ 


-. * Rounded off to nearest 100 psi. 
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Table 9 Endurance strength calculations for the stress ratio 
of 1.0 using the staircase method at 2.5 x 10^ cycles 
for AISI 4340 steel R c 35/40 Phase I grooved 
specimens . 


1 Alternating 
stress 
psi 

i 

n. 

i 

Failures 

in, 

i 

.2 

i n. 

l 






60,168 

5 

i 

5 

25 • 

58,707 

4 

2 

8 

32 

57,246 

3 

6 

18 

54 

55,785 

2 

6 

12 

24 

54,324 

1 

2 

2 

2 

52,863 

0 

1 

0 

0 



N = 18 

A = 45 

B = 137 


d = stress increment = 1,461 psi 

X = lowest stress level = 52,863 psi 
_0 

X = mean (estimate) 

X = X Q + d[A/N " 1/2] = 52,863 ♦ 1,461 [11 - I ] 

X = 55,785 psi a 55,800* psi 
s = standard deviation (estimate) 

s = 1.620 d[(NB-A 2 ) /N 2 + 0.029] « 1.620 (1,461) [ 18(137) - (45) 2 + 0>029 j 
s = 3,290 psi = 3,300* psi (18) 2 


* Rounded off to nearest 100 psi. 





Table 10 Summary of Phase I endurance strength results at 2.5 x 10° cycles 

of life for AISI 4340 steel R c 35/40 grooved specimens for 
distributional Goodman diagram. 
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Table 11 


Diameter and surface hardness at test section of 
AISI 4340 steel Phase II ungrooved specimens 


107 


Specimen 

no. 

Diameter 

in. 

Hardness 
Rockwell C 

1 

0.4977 

37.6 

2 

0.4980 

36.3 

3 

0.4981 

37.0 

4 

0.4979 

36.5 ! 

5 

0.4980 

37.5 

6 

0.4967 

38.3 j 

7 

0.4986 

36.2 

8 

0.4981 

37.0 

9 

0.4979 

37.3 

10 

0.4985 

36.5 

11 

0.4981 

36.5 

12 

0.4983 

36.8 

13 

0.4974 

36.3 

14 

0.4982 

36.3 

15 

0.4983 

37.0 

16 

0.4980 

36.5 

17 

0.4981 

37.3 

18 

0.4982 

36.3 

19 

0.4993 

37.4 

20 

0.4988 

37.6 

21 

0.4986 

39.2 

22 

0 . 4993 

37.5 

23 

0.4980 

36.8 

24 

0.4992 

36.7 

25 

0.4996 

38.3 

26 

0.4979 

37.4 

27 

0.5000 

37.6 

28 

0.5000 

37.0 

29 

0.4981 

36.5 

30 

0.4978 

37.0 

31 

0.4978 

37.6 

32 

0.4983 

37.5 

33 

0.4977 

37.3 

34 

0.4984 

| 37.0 

35 

0.4978 

35.0 


Diameter 


Hardness 


Mean = 0.4983 in. 

Standard deviation = 0.0006 in. 


Mean = 37.05 R c 

Standard deviation = 0.75 R c 




Table 12 


Diameter at the base of groove, groove radius 
and surface hardness of AISI 4340 steel 
Phase II grooved specimens. 
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i "" 

Specimen 

no. 

Diameter 

in. 

Radius 

in. 

Surface hardness 
Rockwell C 

i 

0.4928 

0.35 

38.1 

2 

0.4926 

0.35 

37.6 

3 

0.4909 

0.35 

37.6 

4 

0.4922 

0.35 

38.1 

5 

0.4953 

0.35 

37.3 

6 

0.4966 

0.35 

37.5 

7 

0.4966 

0.35 

36.3 

8 

0.4981 

0.35 

37.6 

9 

0.4936 

0.35 

38.0 

10 

0.4937 

0.35 

37.4 

11 

0.4972 

0.35 

37.5 

12 

0.4938 

0.35 

36.6 

13 

0.4947 

0.35 

37.6 

14 

0.4935 

0.35 

37.5 

15 

0.4973 

0.35 

36.5 

16 

0.4943 

0.35 

37.1 

17 

0.4907 

0.35 

37.2 

18 

0.4948 

0.35 

37.4 

19 

0.4952 

0.35 

37.5 

20 

0.4941 

0.35 

37.0 

21 

0.4969 

0.35 

38.3 

22 

0.4940 

0.35 

36.5 

23 

0.4938 

0.35 

37.5 

24 

0.4930 

0.35 

37.3 

25 

0.4934 

0.35 

37.4 

26 

0.4931 

0.35 

38.0 

27 

0.4937 

0.35 

36.3 

28 

0.4950 

0.35 

37.3 

29 

0.4934 

0.35 

37.2 

30 

0.4924 

0.35 

37.0 

31 

0.4932 

0.35 

36.5 

32 

0.4927 

0.35 

38.1 

33 

0.4939 

0.35 

36.9 

34 

0.4945 

0.35 

37.8 

35 

0.4933 

0.35 

37.0 


Diameter 


Hardness 


Mean = 0.4939 in. 

Standard deviation = 0.0021 in. 


Mean = 37.29 R c 

Standard deviation * 0.53 R c 




Table 13 
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Hardness measurements at Section A of Phase II 
grooved specimens in R £ units see Fig. 17 for locations. 
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Table 14 Hardness measurements for Phase II specimens at 
Section A shown in Fig. 18 in R c units 


Ungrooved Specimens 

Grooved Specimens 

Specimen 

no. 

R c 

units 


X 

J 

Specimen 

no. 

R 

c 

units 


X 

X 

31 

C = 36.5 

36.5 

37.0 

36.7 


33 

C = 36.5 

37.5 

37.0 

37.0 



I = 36.0 

36.0 

35.0 

35.7 



I * 36.0 

38.0 

37.0 

37.0 



0 = 34.5 

35.0 

36.5 

35.3 

35.9 


0 = 36.0 

36.0 

36.0 

36.0 

36.7 

33 

C = 37.5 

38.0 

38.0 

37.8 


34 

C = 35.0 

36.5 

36.0 

36.0 



I = 37.0 

38.0 

37.5 

37.5 



I = 36.0 

35.0 

36.0 

35.7 



0 = 36.5 

36.0 

35.0 

36.0 

37.1 


0 = 37.0 

35.0 

35.0 

35.7 

35.8 

34 

C = 33.5 

35.5 

34.5 

34.5 


35 

C = 37.0 

37.5 

37.0 

37.2 



I = 37.0 

36.0 

37.0 

36.5 


! 

I = 35.5 

37.0 

36.0 

36.2 



0 = 35.5 

35.0 

36.0 

35.5 

35.5 


0 = 35.0 

35.5 

35.0 

35.2 

36.2 
































Table 15 . Tensile test data for Phase II AISI 4340 steel R 35/40 ungrooved specimens. 
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Table 15 Tensile test data for Phase II AISI 4340 steel R c 35/40 ungrooved specimens (Cont'd.) 
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Table 16 


Elongation test data for Phase II AISI 4340 
steel R c 35/40 ungrooved specimens . 
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Specimen 

no. 

Elongation 

2-in. 

gage length, in. 

Elongation 

% 

1 

2.365 

18.25 

2 

2.362 

18.10 

3 

2.365 

18.25 

4 

2.380 

19.00 

5 


17.80 

6 

2.357 

17.85 

7 

2.349 

17.45 

8 

2.368 

18.40 

9 

2.362 

18.10 

10 

2.357 

17.85 

11 

2.357 

17.85 

12 

2.360 

18.00 

13 

2.361 

18.05 

14 

2.361 

18.05 

15 

2.359 

17.45 

16 

2.353 

17.65 

17 

2.354 

17.70 

18 

2.347 

17.35 

19 

2.360 

18.00 

20 

2.361 

18.05 

21 

2.360 

18.00 

22 

2.347 

17.35 

23 

2.340 

17.00 

24 

2.360 

18.00 

25 

2.353 

1 17.65 

26 

2.361 

18.05 

27 

2.358 

17.90 

28 

2.365 

18.25 

29 

2.365 

18.25 

30 

2.361 

18.05 

31 

2.356 

17.80 

32 

2.356 

17.80 

33 

2.356 

17.80 

34 

2.361 

18.05 









Table 17 


Tensile test data for AISI 4340 steel R £ 35/40 
grooved specimens . 


Specimen 

number 

Original 

diameter 

in. 

Original 
cross 
section 
area, in. 

Ultimate 

load 

lb 

Ultimate 

stress 

psi 

1 

0.4928 

0.1906 

52,500 

275,400 

2 

0.4926 

0.1906 

51,950 

272,600 

3 

0.4909 

0.1892 

51,750 

273,500 

4 

0.4922 

0.1903 

51,900 

272,700 

5 

0.4953 

0.1927 

52,200 

270,900 

6 

0.4966 

0.1937 

51,400 

265,400 

7 

0.4966 

0.1937 

51,300 

264,800 

8 

0.4981 

0.1949 

52,200 

267,800 

9 

0.4936 

0.1913 

52,000 

271,800 

10 

0.4937 

0.1914 

51,500 

269,000 

11 

0.4972 

0.1942 

52,300 

269,300 

12 

0.4938 

0.1915 

51,600 

269,500 

13 

0.4947 

0.1922 

52,000 

270,600 

14 

0.4935 

0.1912 

51,600 

269,900 

15 

0.4973 

0.1942 

51,900 

267,300 

16 

0.4943 

0.1919 

51,450 

268,100 

17 

0.4907 

0.1891 

51,300 

271,300 

18 

0.4948 

0.1923 

51,600 

268,300 

19 

0.4952 

0.1926 

52,100 

270,500 

20 

0.4941 

0.1917 

51,600 

269,200 

21 

0.4969 

0.1939 

52,500 

270,800 ! 

22 

0.4940 

0.1916 

51,250 

267,500 

23 

0.4938 

0.1914 

51,150 

267,200 

24 

0.4930 

0.1909 

51,050 

267,400 

25 

0.4934 

0.1912 

51,800 

270,900 

26 

0.4931 

0.1909 

51,650 

270,600 

27 

0.4937 

0.1914 

50,900 

265,900 

28 

0.4950 

0.1924 

51,550 

267,900 

29 

0.4934 

0.1912 

51,500 

269,400 

30 

0.4924 

0.1905 

51,000 

267,700 

31 

0.4923 

0.1904 

49,600 

260,500 

32 

0.4927 

0.1907 

1 

51,250 

268,700 


Summary of static strength results 
grooved specimens. 

Ultimate strength 


Mean 


= 269,000 psi 


Standard deviation = 2,800 psi 























Table 19 


Calibration coefficients and speed for each research 
machine and for Mode 5 operation*. 


Machine 

no. 

K 

BGR 

K 

GR-TH 

K 

T 

K 

T/B 

K 

B/T 

Machine 

speed 

rpm 


1.3081 

0.0149 

0.8753 


0.0477 

1779 

2 

1.3081 

0.0145 

0. 8950 

0.0463 

0.0555 

1775 

3 

1.3081 

0.0170 

0.8748 

0.0542 

0.0645 

1778 


*Mode 5 operation refers to the computer program code which identifies 
the proper calibration coefficients to be used with each set of machine 
data obtained during a specific calendar period. Mode 5 calibration 
coefficients should be used for all data obtained after June 1, 1971 until 
the next calibration. 
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d = stress increment = 1,900 p S i 

X = lowest stress level = 77,400 psi 
o 

X = mean (estimate) 

X » X Q + d[A/N - 1/2] =, 77,400 + 1,900 [1£ _ I ] 

X= 80,725 psi 2 80,700* psi 13 2 

s = standard deviation (estimate) 

S = 1.620 d [ (NB-A 2 ) + 0.029] * 1.620 (1,900) f 13(77) - (29) 2 n n?Q ] 

S = 3,040 psi 5 3,000* psi ^ .2 


Rounded off to nearest 100 psi. 
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Table 21 Endurance strength distribution parameters calculations 
for 0..25QQ in. Phase II research AISI 4340 steel, R c 
35/40 ungrooved specimens for stress ratio of <*>. 


Alternating 

stress 

psi 

i 

n. 

i 

Failures 

in. 

i 

.2 

1 n i 






83,100 

3 

3 

9 

27 

81,200 

2 

9 

. 18 

36 

79,300 

1 

4 

4 

4 

77,400 

0 

0 

0 

0 



N = 16 

A = 31 

B = 67 


d = stress increment =1,900 psi 
X = lowest stress level = 77,400 psi 
X = mean (estimate) 

X = X Q + d [A/N - 1/2] = 77,400 + 1,900 [li - i ] 

X = 80,230 psi ^ 80,200* psi 16 2 

s = standard deviation' (estimate) 

S = 1.620 d[(NB-A 2 )/N 2 + 0.029] = 1.620 (1,900) £ 16(67) - (31) 2 + Q Q2g 

s = 1,425 psi 2 1,400* psi (16) 2 


* Rounded off to nearest 100 psi. 
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Table 23 Normal cycles-to-failure distribution parameters for Phase II results. 
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O P< 
o 


D critical at 0.05 level of significance and sample size of 35 is 0.224. 
Accept/reject critical value for chi-squared test at 0.05 level of significance: 
Degrees of freedom Critical value 

1 3.841 

2 5.991 
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B = 88 


d = stress increment = 2,053 psi 

X = lowest stress level = 28,762 psi 
o 

X = mean (estimate) 

X = X q + d[A/N + 1/2] =. 28,762 + 2,053 [ 

X = 33,689 psi 2 33,700* psi 
s = standard deviation (estimate) 


38 1 i 

— + — J 
20 2 


S = 1.620 d [ (NB-A j /hi + 0.029]= 1.620 (2,053) [ 20(88) - (38) 0>02 9] 
s = 2,724 psi s 2,700* psi (20) 2 


* Rounded off to nearest 100 psi. 













Table 26 Endurance strength calculations for Phase II grooved 
specimens at stress ratio of 1.06. 


Alternating 
s. tress 



Failure' 



35,202 

34,109 

33,016 

31,923 


30,830 

29,737 



d * stress increment = 1,093 psi average 
X Q = lo\;est stress level = 29,737 psi 
X = mean (.estimate) 

X = X q + d[A/N - 1/2] = 29,737 + 1,093 [51 0.5 ] 

X = 31,978 psi 2 32,000 * psi 20 


s - standard deviation' (estimate) 

' 1.620 d[(NB-A 2 )/N 2 + 0.029] = 1.620 (1,093) [ 
: 2,796 psi 2 2,800 *psi 


20(161 


Rounded off to nearest 100 psi. 
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Table 27 ’ Endurance strength calculations for Phase II 
grooved specimens at stress ratio of 0.40. 


Alternating 

stress 

psi 


33.000 

32.000 

31.000 

30.000 

29.000 

28.000 




d = stress increment = 1,000 psi 

X = lowest stress level = 28,000 psi 
__o 

X = mean (estimate) 

• X = X + d[A/N + 1/2] = 28,000 + l,000[fi i ] 

— ° 22 2 
X = 34,000 psi 

s = standard deviation' (estimate) 

S = 1.620 d[ (NB-A 2 ) /N 2 + 0.029] = 1.620 (1,000) r 22(115) - (41) 2 + Q Q2g ] 

5 = 2,900* * psi (22)2 - 1 


* Rounded off to nearest 100 psi. 















Table 28 Endurance strength calculations for Phase II 
grooved specimens at stress ratio of 0.25. 


Alternating 

stress 

psi 



24,000 


23,000 


d = stress increment = 1,000 psi 

X = lowest stress level = 23,000 psi 
o 

X = mean (estimate) 

X = X q + d [A/N - 1/2] = 23,000 + 1,000[ i£ _ l ] 

X= 24,136 psi 2 24,100 * psi 22 2 

£ = standard deviation' (estimate) 

S = 1.620 d[(NB-A 2 )/N 2 + 0.029] = 1,620 (1,000) \ 22(84) - (36) 2 + 0 029 ] 
s- 1,732 psi 5 1,700* psi L (22) 2 + ’ J 


* Rounded off to nearest 100 psi. 
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d = stress increment =1,000 psi 

X = lowest stress level = 22,000 psi 
__o 

X = mean (estimate) 


X = X q + o [A/'N - 1/2] =. 22,000 + 1,000 [ ££ _ £ ] 

X = 22,900^ s i 22 2 

s = standard deviation' (estimate) 

S = 1.620 d((NB-A 2 )/M 2 + 0.029] = 1.620 (1,000) f 22(53) - (31) 2 + Q Q2g ] 
S = 700* psi ^ (22) 2 


* Rounded off to nearest 100 psi. 
















Table 30 Summary of Phase II endurance strength results at 2.5 x 10° cycles of 

life for AISI 4340 steel R c 35/40 grooved specimens for distributional 
Goodman diagram. 
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APPENDICES 


APPENDIX A 


PDP PROGRAM TO CALCULATE VISICORDER DIVISIONS FOR 130 

DESIRED LEVELS OF BENDING AND TORQUE 

C-8K MODV 11-219 

01.09 T "MACHINE NO. 3", !, !, ! 

01.10 T "STRESS RATIO = INFINITY", ! 

01.20 A "STRESS LEVEL = ", SL, ! 

01.30 A "BENDING CAL. RESIST. =", RB, ! 

01.50 S K1 = 1.3081 

01.60 S K2 = .0170 

01.70 S K3 * .8748 

01.80 S K4 = .0542 

01.90 S K5 = .0645 

02.06 S DS = .491 

02.07 S DO = 2.0 

02.08 S DI = 1.3125 

02.20 S NB = 50 

02.30 S NA = 4 

02.40 S GB = 3.23 

02.50 S A = 30000000 

02.60 S BR = 190 

02.70 S NT = 45 

02.80 S GT = 2.06 

02.90 S TR + 120 

03.10 S TO = 0 

03.20 S SO = K2*K1*SL 

03.30 S SP = S0+K4*T0 

03.40 S VB = (NB*NA*GB*RB*SP)/ (A*BR) 

03.80 T %7.03, "VIS. DIV. - BENDING", VB, ! 

03.90 T %8. 01 , "BENDING CALIB.", RB," FOR 50 DIV.,", ! 

05.10 S EB = . 03*VB 

05.20 S ET + . 03*VT 

05.30 S SB + SL/VB 

05.50 T %5. 02 "DIV. OF BENDING FOR 3% ERROR", EB, » 

05.70 T %8. 01 "BENDING SENSITIVITY", SB," PSI/DIV.", ! 
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APPENDIX B 

PROGRAM STRESS (Updated) 

The purpose of this program is to convert each Visicorder record 
into normal stress, shear stress and stress ratio. The program 
calculates the mean and standard deviation of these two stresses and 
of the stress ratio of each group of cycles-to-failure data. It 
calculates the cycles-to-failure from time-to-failure data. The 
program will also accept cycles-to-failure data. 

This program remains basically as described previously. 

Added are the latest calibration coefficients and a routine to 
convert shear stresses to normal mean stresses. 

The program distinguishes between cycles-to-failure and time-to- 
failure through the use of a code number. The code also tells the 
program whether or not the group of data are endurance test data. This 
discrimination is necessary because the program will not calculate the 
endurance strength distribution parameters. The discrimination code 
is fed in as data, and is as follows: 

0 if the failure data is in time to failure. 

1 if the data is endurance test data. 

2 if the failure data is in cycles to failure. 
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The program will accept as many sets of data as desired and the 
groups may be mixed; i.e., endurance test, group with cycles-to-failure 
data, and group with time to failure data. A group is all the data 
for one stress level. A group of data consists of the following: The 

first card contains in this order, the number of specimens in the level, 
the mode of the run and the code. The mode is dependent upon the date 
the run was made and identifies the calibration constants for each 
machine to be used in the computations. 

The fields on the data card are as follows: 
spaces 1 to 5 - number of specimens 
spaces 6 to 10 - mode 
spaces 11 to 15 - code. 

The number of specimens, mode and code are fixed point numbers and have 
no decimals but the numbers must be placed to the right in each field. 

The next sequence of cards reads in the cycles-to-failure or time 
to failure in hours, minutes and seconds. If the data is in time to 
failure there are ten groups on a card, so the number of cards required 
will depend upon how many specimens are in the stress level. The format 
across the card is as follows: 


spaces 

i 

and 

2 

- blank 

spaces 

3 

and 

4 

- hours 

spaces 

5 

and 

6 

- minutes 

spaces 

7 

and 

8 

- seconds 

spaces 

9 

and 

10 

- blank > 
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and the sequence continues in this manner. If the failure data is in 
cycles -to-fai lure the format is 8 fields of 10 spaces each and the 
decimals appear in the last space of each field; i.e. spaces 10, 20, 

30, etc. If the group of data is for an endurance test there is no 
failure data and these cards are left out. The program will automatically 
' handle the data if the proper code number is put on the first card. 
Following the cycles-to-failure cards are the cards containing the 
information for each specimen in the stress level . The information 
must be placed on each card as follows: 
spaces 1 to 5 - test number 

spaces 6 to 10 - specimen number 

spaces 11 to 15 - machine number 

spaces 16 to 20 - with a decimal in space 20 - pan weight 

spaces 21 to 30 with a decimal in space 28 - bending calibration 

resistance 

spaces 31 to 40 with a decimal in space 38 - number of bending 

calibration divisions 

spaces 41 to 50 with a 

decimal in space 48 - number of divisions of bending 
spaces 51 to 60 with a 

decimal in space 58 - torque calibration resistance 

spaces 61 to 70 with a 

decimal in space 68 - number of torque calibration divisions 



spaces 71 to 80 with a 

decimal in space 78 - number of divisions of torque 

The test number, specimen number, and machine number are fixed 
point numbers and must be placed to the right in each field. There 
is one card for each test specimen and the cards must be placed in 
the same order as the failure data is placed on the cards preceeding 
these cards. For data at stress ratio of ® there will be no torque 
stress data. In this case these fields can be left blank. The 
computer reads blanks on data cards as zeros. 

This makes up one group of data at a given stress level and 
ratio. As many groups may be run as desired by simply placing the 
groups one behind the other. 

A list of important variables and symbols in the program STRESS 
using Fortran language follows: 



NCARDS 

MODE 

NCODE 


XHOURS(I) 
XMIN(I) 
SECS (I) 

TOTCY(I) 

NOTEST 

NOSPEC 

MACHNO 

PANWT 

RCALB 

ENCLAB 

ENVISB 

RCALT 

ENCALT 

ENVIST 

ENA 


List of Definitions for Program to Find 
Stress Levels and Ratios (PROGRAM STRESS) 

number of specimens tested at given level, 
number of mode depending on date of test. 

0 if failure data is in times to failure. 

1 if data is from an endurance level . 

2 if data is in cycles-to-failure . 

times to failure in hours, minutes and seconds. 

cycles-to-failure . 
test number, 
specimen number, 
machine number. 

amount of weight on loading arm. 

calibration resistance used in bending channel. 

number of visicorder divisions used when calibrating 
bending channel . 

number of divisions during actual test. 

calibration resistance used in torque channel . 

number of visicorder divisions used when calibrating 
torque channel . 

number of divisions during actual test, 
number of active arms in strain gage bridge. 



RGAGEB 

= 

resistance of bending strain gage. 

RGAGET 


resistance of torque strain gages. 

GB 

= 

bending gage 

factor. 

GT 

= 

torque gage 

factor. 

CBGR 

= 

calibration 

constant Kgg R . 

CGRTH 

- 

calibration 

constant K GR _ TH - 

CT 

= 

calibration 

constant K^. 

CTB 

= 

calibration 

constant K^g. 

RPM 

= 

revolutions 

per minute of machine. 

SOUTH 

S 

output normal stress corrected for interaction. 

TAUTH 

= 

output shear 

stress corrected for interaction. 

STRGR(I) 

= 

normal stress in specimen groove. 

TAUGER(I) 


shear stress 

in specimen groove. 

SOUTHP 

= 

output stress not corrected for interaction. 

TAUTHP 

= 

output stress not corrected for interaction. 


The program STRESS listing in Fortran language follows: 
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PROGRAM STRESS ( INPUT, OUTPUT, TAPE1 =INPUT) 

DIMENSION STRGRC50), TAUGR150), R (50) ,XHOURS (50 ) ,XMIN(50 ) » 
1SECS(50) ,TOTCY<50) 


C READ IN THE NUM9ER OF CAROS IN THE STRESS LEVEL ANO THE MODE 

C OF OPERATION OF THE MACHINE 

C NCOOE = 1 FOR ENDURANCE LEVEL, NCOOE = 0 IF FAILURES IN TIMES TO FAILUft 

C NCODE = 2 IF FAILURES IN CYCLES TO FAILURE. 

35 READ 10C, NCAPQS, MODE, NCOOE 

100 FORMAT (315) 

IF ( EOF , 1) 201, 171 

171 PRINT 91 

91 FORMAT (1H1//) 

170 IF (NCODE. EO.l) GO TO 175 

C FORMAT FOR OUTPUT HEADINGS INCLUOING TIMES TO FAILURE. 

PRINT 310, MODE 

310 FORMAT (62X,11HTEST MODE = , I 2/ ) 

PRINT 61 

61 FORMAT (2X,34HTEST SPEC. MACH. PAN CYCLE S , 7X , 4HRC AL , 7 X 


1 , 4 HNGAL ,7 X , 4HN VIS , 6X , 4HRCA L , 6X , 4HNCAL , 6X , 4HNV IS , 5X , 7H9EN0ING »4X » 1 
26HSHEAR STRESS/ 3X , 35HNO . NO. NO. WT . TO FAILURE, 3X, 

37H3ENQING, 4X , 7 HBENDING , 4 X , 7HBE NO I NG , 4X , 6HT OR DUE , 4X, 6HT0RGUE , 4X , 

46H TO R0UE,5X,6HSTRESS,4X,6HSTR£SS,sX,5H RATIO//) 

IP (NCODE. £0.2) GO TO 320 


C ROUTINE TO READ IN TIMES TO FAILURE. 

READ 401, (XHOURS (I) ,XMIN(I> , SEC 3(1) , 1=1, NCAfiOS) 
401 FORMAT ( 1 0 ( 2X , 3F2 . 0 ) ) 

GO TO 30 0 

O READ IN CYCLES TO FAILURE. 

320 READ 402, (TOTCY(I), I=1,NCAPDS> 

402 FORMAT ( 8F1Q.0) 

GO TO 300 

C FORMAT FOR OUTPUT HEADINGS WITHOUT TIMES TO FAILURE. 

175 PRINT 90, MODE 

90 FORMAT ( 5 0 X , 1 <*HENQUP ANCE TEST , 10X , 11HTEST MODE =,I2/> 
172 PRINT 6 C 


60 FORMAT (3X,4HTEST,3X,8HSP£CIMEN,3X,7HMACHINE , 5X , 3HPAN , 7X ,4 HRCA L , 7X 
1 , 4HNCA L , 7 X , 4HNV IS , 6X , 4HRC AL ,6X , 4HNCAL , 6X , 4HNV IS , 5 X , 7H9ENDI NG ,4X , 


25 HSHEA R , 4X ,6HSTRESS/4X, 3HNO . , 6X, 3HN0. , 7 X , 3HN0. , 5 X , 6HW£ IGHT , 4X , 
37H8ENOING , 4X , 7H9EN0ING , 4X , 7Hfi£ND ING, 4X , 6HT ORQUE , 4 X , 6HTOR QUE , 4X , 
46HT0RQUE , 5 X , 6HSTRESS , 4X , 6HS T RSSS , , 5HR AT 1 0/ ) 

C ROUTINE TO CALCULATE POLAR MOMENTS OF INERTIA. 

C SPOIA = SPECIMEN DIA. TOO = TOOLHOLOER O. O. TIO = TOOLHOLOSR I. O 


300 SP0IA=.491G 
TO D=2. 0 
TI 0=1. 3125 

SPJC=( 3. 14159*SPDIA**3) /16. 0 

THJC = (3.1«.159*(T0D**4-TID**4))/(16.*T0D> 

CONST=THJC/SPJC 
DO 120 I = 1, NCARDS 

C READ IN THE TEST NO., SPECIMEN NO., MACHINE NO., PANWEIGHT, AND THE 



O O -H o O O 
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CALIBRATION RESISTANCE', VISICORDER CALIBRATION DISTANCE AND VIS ICO, ROE R 
OUTPUT DIVISIONS FOR BENDING ANO TORQUE 
RO READ 10, NOTEST, NOSPEC ,H ACHNO , PANHT , RCAL6 , ENCAL B ,ENV IS B , RC A L T , 
1ENCALT.ENVIST 

0 FORMAT C 315 ,F5. 0 ,6F10. 2) 

CAROS = NC ARDS 

OEFINE THE ELASTIC MODULUS, NO. OF ACTIVE ARMS OF THE BRIDGES, THE 
RESISTANCES OF THE BENDING ANO TORQUE GAUGES, AND THE BENUING ANO 
TORQUE GAUGE FACTORS 
E s 30000000* 

ENA=4. 

RGAGEB=19Q . 

RGAGET=120. 

G8- ?• 2 3 

GT = 2. 06 

C SELECTION OF MACHINE ANO MODE 

IFCMACHNO. EO.l. AND. MODE. EQ.l) GO TO 11 
IFCMACHNO. EQ.2.AND.MODE.EQ. 1) GO TO 24 
IF(MACHNO.E0.3.ANO.MODE.E0.1) GO TO 31 
IF (MAC HNO. EQ.l. AND. MOOE.EO.2) GO TO 11 
IFCMACHNO. EQ. 2. ANO. MODE. EQ. 2) GO TO 24 
IFCMACHNO. EQ. 3. AND. MODE. EQ. 2) GO TO 31 
IFCMACHNO. EQ.l. AND. MODE. EQ. 3) GO TO 11 
IFCMACHNO. EQ. 2. ANO. MODS. EQ. 31 GO TO 24 
IFCMACHNO. EQ. 3. ANO. MODE. EQ. 3) GO TO 33 
IFCMACHNO. EQ.l. ANO. MODE. EO. 4) GO TO 11 
IFCMACHNO. EQ. 2. ANO.MODE. EO. 4) GO TO 24 
IFCMACHNO. £0.3. ANO. MODE. EQ. 4) GO TO 34 


IFCMACHNO . EQ . 1 

.AND. 

MODE .EO. 

5) 

GO TO 

15 

IFCMACHNO . EO. 2 

. ANO. 

MODE .EQ. 

5) 

GO TO 

25 

IFCMACHNO .EO. 3 

. ANO. 

MODE . EQ. 

5) 

GO TO 

53 

CALIBRATION param 

ETERS 

FOR GIVEN 

MOD 

£ ANO 

MACHINE 


11 CBGR = 1.0123 
CGRTH = .0208 
CT = .3752 
CT9 = -.0459 
C3T = .029 
RPM=1786. 

GO TO 50 
15 CBGR =1.3031 
CGRTH=0. 0149 
CT =0.87g3 
CT8=-0.0462 
CBT=0. 0493 
PPM = 1779 
GO TO 50 

24 CBGR = 1.0123 
CGRTH = .0183 
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CT = .8201 
CTP = 0.0344 
C3T = .0422 
RPM=1784. 

GO TO 50 

25 C3GR- 1.3081 
CGR7H= 0.0145 
CT = 0 • 895 0 

' CTR=0.0463 

C3T=0. 0544 
R9M = 1775 

, GO TO 5 G 

,'31 CRGR = 1.0945 
CGRTH = .0211 
CT = .933 
CT 8 = .G 
CRT = -.0149 
R°M = 17 8 C . 

GO TO 50 

33 CRGR = 1.0945 
CGRTH = .0211 
CT - .7721 
CTB = .0 

CRT = -.0127 
RPM= 17 3 0 . 

. GO TO 50 

34 CBGR = 1.0123 
CGPTH = .3197 
CT = .7721 
CTR = .0 

C3T = -.0127 
RPM-1780. 

53 CBGR= 1 . 3 03 1 
CGP.THs0.Gl7C 
CT= 3.8748 
CT R = 0. 0542 
CBT = C. 0645 
RPH s 1773 
GO TO 50 

50 IF (ENVIST.EO.0.0) GO TO 160 
GO TO 51 

C CALCULATION OF RENOING STRESS LEVEL FOR INFINITY RATIO 

160 SOUTH=(ENVISR*E*PGAGFR) / (ENCALR*ENA*G8*RCAL8> 

STRGR(I) = SOUTH / (CGRTH * CRGR) 

TAUGR(I) = 0.0 
IF (NCODE.EO.l) GO TO 301 
IF (NCODE.EO.2) GO TO 330 

C CALCULATE CYCLES TO FAILURE FROM TIMES TO FAILURE. 

CY HR = XHOURS( I) ♦60.*RPM 
CYMIN = X M IN ( I ) *RPM 
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CYSEC = (SECS<I>*RPM>/60.U 
TCY = CYHR+CYMIN+CYSEC 

GO TO 335 

330 TCY = TOTCY(I) 

3 35 PRINT 7 2, NOTE ST, NOS°E C , MACH NO , FA N WT , TCY, RC AL 3 , c NCA L8 , E N V IS 8 , RCALT, 
1ENICALT,ENVIST,STRGR( T) ,TAUGR<I> 

7? FORMAT (3I6,P5. 1,2P11. 0,F10, 2,F11.2,F11. 0 ,F1C.2, 

lF10.2, r 11.0,F10.0,4X,6HINRIN./> 

GO TO 120 

301 PRINT 71, NOTEST, NOSPEC , MACHNO , PANWT ,RCALB ,ENCAL 8,c.NVIS3 ,RCALT , 
1ENCALT , ENV 1ST ,STRGR(I) ,TAUGR(I> 

71 FORMAT (4X, I3,&X,I3,9X,H,6X,F5. 1, F12. 0,F10.2,F11.2,FH. C,F10.2, 

1F1G.2,F11. C,FlO.O,3X,fHINFIN. /) 

GO TO 120 

C CALCULATION OF RENDING STRESS, SHEAR STRESS AND STRESS RATIO FOR 

C ALL FINITE RATIOS 

51 SOUTHP=<ENVISR*E»RGAGFB) / ( ENCAL3»ENA*GB*RCALR) 

TAUTHP= < EN VIST *E*RG A GET )/< ENC ALT* SNA *GT* RCALT ) 

SO UT H= SO UT HP -CTR*TAUT HP 

TAUTH=TAUTHP-CRT*SOUTH 

STRGR(I) = SOUTH / (CGPT H * C3GR) 

TAUGR(I) = CT*TAUTH*CONST 
R ( I ) = STRGR(I) / (TAUGR(I) * 1.772) 

IF (NCODE.EO.l) GO TO 302 
IF (NC0CE.EQ.2) GO TO 3UC 

C CALCULATE CYCLES TO FAILURE FROM TIMES TO FAILURE. 

CYHR = XHOURS(I)*60 .*RPM 
CYMIN = XM IN ( I) *RPM 
CYSEC = (SECS(I)*PPM)/6C.Q 
TCY = CYHR+CYMIN+CYSEC 

GO TO 3<*5 

340 TCY = TOTCY(I) 

345 PRINT 7 3, NOTEST, NOS PEC, MACH NO, PANWT, TCY, RC AL 3 , E NCAL 3 , E NV IS 8 , RCALT, 
ltNCALT , ENV 1ST ,STRGR(I) , T AUGR ( I ) ,R(I) 

73 FORMAT ( 3 16 , FA . 1 , 2 FI 1 . 0 , Fl J . 2 , FH . 2 , fu . C , FI 0 . 2 , 
1F10.2,F11.0,F10.0,F9.3/) 

GO TO 120 

302 PRINT 70 , NOTEST, NOSPEC, MACHNO, PANWT , RCALB , ENCA L O , EN VIS 3 , RCA LT , 
1SNCALT ,ENV 1ST »ST p GR(I) ,TAUGR(I),R(I) 

70 FORMAT (4X,I3,6X,T3,9X,I1,oX,F5.1,F12.0,F1C.2,F11.2,F11.G,F10.2, 

1F10.2,F11.0,F10.0,F9.3/) 

120 CONTINUE 

IF (NCODE.EQ.I.ANO.ENVIST.EO.O.O) GO TO 35 
IF (NCODE. EQ.l) GO TO 200 

C CALCULATION OF MEAN AND S T AN DAP 6 Q£ V I A T I O N OF R.ENOING STRESS, SHEAR 

C STRESS, AND STRESS RATIO 

CALL MEAN (STRGR , CARDS ,NCAPDS, XME AN , OEV) 

PRINT 3 

3 FORMAT (1H ) 

PRINT 130, XME AN , OEV 
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130 FORMAT <19X,3iHMEAN BENDING STRESS IN GROOVE =,F1Q.1,5H PSI.//11X, 
139HSTD. DEV. OF BENDING STRESS IN GROOVE =,Fl0.2,5H PSI./) 

132 IF(ENVIST.EQ.O.O) GO TO 35 
IP (NCODE.EQ.l) GO TO 200 

CALL MEAN (TAUGR, CARDS, NCAROS, XMEAN , DEV) 

PRINT 3 

PRINT 140, XMEAN » DEV 

140 FORMAT (20X,30HMEAN TORQUE STRESS IN GROOVE =,F1G.1,5H ^SI./ 12X, 
133HSTD. DEV. OF TORQUE STRESS IN GROOVE =,fiq.2,5H PSI. ) 

SMEAN = XMEANM.732 
SDE V = DEV*i. 732 
PRINT 3 

PRINT 142, SMEAN, SDEV 

1^2 FORMAT (19X , 'MEAN NORMAL STRESSSIN GPOOVE*, F10.1, * PSI.*//11X, 
i*STQ. DEV. OF NORMAL MEAN STRESS IN GR = *, F10.2, * P3I.+/I 

200 s CALL MEAN ( R, CAROS , NCAROS , XMEAN, DEV) 

PRINT 3 

PRINT 150, XMEAN, DEV 

150 FORMAT ( 3 IX , 19HMEAN STRESS RATIO =,F10.5/ 23X,27HSTO, DEV. OF S TR£ 
1SS RATIO =,F10.5> 

37 GO TO 35 

201 STOP 
END 

SUBROUTINE MEAN (X, DATA, NOATA, XMEAN, OLV) 

C SUBROUTINE TO CALCULATE THE MEAN AND STANDARD DEVIATION OF DATA. 

DIMENSION X(NDATA) 

S I G M A = 0.0 
DO 3 1=1, NDATA 
tt SIGMa=SIGMA+ X(I> 

XMEAN = SIGMA/DATA 
TOP2 = 0.0 
DO 9 1=1, NDATA 

9 T OP 2 = T OP 2 + (X(I) - XMEAN)*»2 

OEV =SQPT (TOP2/CDATA - 1.0)) 

RETURN 

.END 
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APPENDIX C 

PROGRAM CYTOFR (Updated) 

This program calculates estimates of the mean and standard 
deviation of the cycles -to-failure data for both the normal and the 
lognormal distributions, and calculates the moment coefficients of 
skewness and kurtosis. It also performs the Chi-squared and the 
Kolmogorov -Smirnov goodness-of-fit tests. 

The program has been updated since it was previously reported. 

The main program was revised to incorporate a sort 
routine to preclude the necessity for manually ordering the data 
inputs. Subroutine CHISQA (Chi-squared test) was modified to provide 
for automatic combining of cells at the tails of the distribution 
when the end cells do not contain at least five failure data points. 
Subroutine GRAPH was added to the program to plot a histogram of the 
cycles-to-failure data based on the cell widths and number of 
failure data points per cell computed by the CHISQA subroutine. The 
theoretical distribution curve represented by the parameters 
estimated by the main program is sketched and superimposed over the 
histogram. The plotting of the histogram and of the distribution is 
done by the Cal-Comp plotter from an output tape generated by the 
computer. 

The data deck for operating program CYTOFR is in three logical 
sections per problem. The first section consists of two cards. Card 
one consists of three titles, while card two specifies the parameters 
necessary for the statistical calculations in CYTOFR. 
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The second, logical section is a variable number of cards, each 
specifying up to eight data points for analysis. 

Section three provides parameters for plotting the normal and 
lognormal distributions. The first card is a parameter list for the 
normal distribution plot, followed by exactly six cards of footnotes. 
Next is a parameter list card for the lognormal distribution. 

Multiple problems may be executed by stacking complete data sets 
behind each other. A list of important variables and symbols in 
program CYTOFR using Fortran language follows: 



List of Definitions for Program to fit Normal 
and Log-Normal Distributions to Cycles - 
to-Failure Data (PROGRAM CYTOFR) 


Main Program: 

NDATA = DATA = number of observations. 

STRLV * stress level in psi. 

AKURCY = accuracy to which cycles-to-failure data are known. 

RATIO = stress ratio 

X(I) = cycles-to-failure data 

CUMFRQ(I) = cumulative frequency of each X(I) ; ie, number of 
X's less than or equal to X(I) . 

PCAREA(I) = CUMFRQ ( I ) /NDATA 


Subroutine to calculate the mean and standard deviation of the 
cycles-to-failure data (SUBROUTINE MEAN) 

SIGMA = sum of the X(I)'s 

XMEAN = average of the X (I) 's 


TOP 2 = £ (X(I)-XMEAN) 2 

i=i 

DEV = standard deviation of the X(I)*s 


Function subroutine to find the area under the normal curve 
(FUNCTION PROB(X) ) . 

X = abscissa value for which corresponding area 

is desired. 


PROB 


desired area. 
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Subroutine for Chi-square goodness -of- fit test (SUB-ROUTINE CHISQA) . 
K = number cells. 

XMAX = largest value of cycles-to-failure. 

XMIN = smallest value of cycles-to-failure. 

CSV = cell starting value. 

CEV = cell end value. 

CLB = cell lower bound. 

CUB = cell upper bound. 

FREQ(J) = number of observations in J^ 1 cell. 

REQAREA(J)= expected value of cell. 

CHISQR = total Chi-square value. 

U(I) = Chi-square value of I** 1 cell. 

Subroutine for Kolraogorov-Smirnov test (SUBROUTINE DTEST) . 

Z(I) = abscissa value on standard normal curve for a 

given X(I) . 

ARUNCN = area under standard normal curve from - to Z(I). 
DSTAT(I) = absolute difference between the data cumulative 

frequency and the hypothesized cumulative frequency. 
XMEAN = average of the X(I)'s. 

DEV = standard deviation of the X(I)'s 

PROB(T) = area under the standard normal curve from -T to 


+T. 
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Subroutine to calculate the moment coefficients of skewness and 
kurtosis (SUBROUTINE ALPHA) . 

ALPHA3 = moment coefficient of skewness. 

ALPHA4 = moment coefficient of skewness, 

n 2 

VAR = Z (X(I)-X) 

i=i 

n . 

T0P3 = E (X(I)-X) 3 

i=i — 

SKEW = third moment of the data. 

STDEV = biased estimator for standard deviation. 

n . 

TOP 4 = Z (X(I)-X) 4 

i=i — 

TKURT = fourth moment of the data. 
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Card 

1 


2 


3 to 
n + 1 


DATA DECK STRUCTURE 


Columns 

1-20 


21-40 


Description 

Twenty character descriptive title 
to appear at the top of each printed 
output page. 

Twenty character descriptive title of 
the input data will appear on printed 
output, as well as the X-axis label for 
both slots. 


41-50 Unit of data measurement (cycles, inches 

etc. ) . 

1-10 Number of data points on following 

card(s) . Must have a decimal point. 


11-20 

21-30 


Stress level. Must have a decimal point 
Stress ratio. Must have a decimal point 


31-40 

.** 1 - 10 , 11 - 20 , ..., 
71-80 

1-10 


Accuracy. Must have a decimal point. 

Data points, punched eight per card 
until all points are exhausted. Must 
have decimal points. 

The letters "NORMAL" followed by four 
blank spaces. 


11-20 Length of the X-axis in inches. If 

zero or blank, 6.0 is assumed. Decimal 
point is necessary. 

21-30 Length of the X-axis in inches. If zero 

or blank, 5.0 is assumed; if greater 
than 5.0, 5.0 is assumed. Decimal point 
is necessary. 


**The notation 3rd to n is intended to mean from the third card to the nth 
card. In the case of 34 data points, the second section of data would 
stretch from the 3rd to 7th card . 
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DATA DECK STRUCTURE (Continued) 


Card Columns 

31-40 

41-50 

51-60 

61-70 


71-80 


n + 2 to 1-50 

n + 7 

n + 8 1-10 


Description 

Minimum X-axis value. The plotting 
program may find it necessary to alter 
this value slightly. If absent or 
zero, a resonable maximum will be 
assumed. The decimal point is necessary. 

Minimum Y-axis value. This should be 
blank or 0. 

Maximum X-axis value. Follow same 
rules as minimum X-axis value. 

Maximum Y-axis value. If absent or zero, 
the plotting program will search for the 
smallest even number greater than Cor 
equal to) the height of the tallest 
histogram block, an automatic adjustment 
will be made. The decimal point is 
necessary. 


Height of lettering on graph, if 0 < 
height < .15. Otherwise this parameter 
will be set to equal .15. 

Footnotes punched as they will appear 
on the normal graph. 

The letters "LOG-NORMAL." Otherwise 
the same as card n + 1 . 


The program for CYTOFR listing in Fortran language follows: 
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PROGRAM CYTOFR (INPUT ♦OUT PUT, TAPE 1= INPUT , TAPE?) 

0 PROGRAM TO FIT NORMAL AND LOG-NORMAL CURVE TO DATA ANO CHECK 

C GOODNESS OF FIT. 

COMMON CUMFRQC100) ,NDATA,X (100) , DE V , XME AN , CLB ( 9 > , CUR ( 9) . FRE 0 ( 9) , K , 

1 CHISQR,TITLE(2> , SUBTIT L (2) ,CS V ( 9 ) , CEV ( 9 ) , PC ARE A ( 10 0 > ,DSTAT(10Q»* 

2 AREA (9) , REQAREA (9) , EXFRE 0 (9 ) ,U ( 9) , Z ( 10 0 ) ,NX ( 1 0 0 > , DATA , AKIJRCY , 

3 XMAX,XMIN,PSI,CD,D,PATIO,XLENGTH,YLENGTH,YMAX, YMIM,XMA,XMI , 

4 HL£TT£R,COM(3) , W, ALPHA3 , A LPHA4 , FOOT (30) ,IT,SKS,UNIT 
INTEGER TITLE, SUBTITL, FOOT 

EXTERNAL °ROB 
C 

C INITIALIZE PLOTTER 

C 

CALL INITIAL (0,2, 0.3,0) 

71C PRINT 1 

C NDATA=DATA-NUMBER OF OBSERVATIONS 

C STRLV = STRESS LEVEL IN PSI. 

C X=NUMRER OF CYCLES TO FAILURE 

READ 50C, TITLE, SUBTITL, UNIT 

500 FORMAT ( 8A1 0) 

IF (EOF (in 56,5 

5 READ 501 , DATA ,STRLEV, RATIO, AKUFCY 

501 FORMAT (BF1C.0) 

NDAT A=OA T4 

READ 501, (X(I) ,1=1, NO AT A) 

C 

C SORT X(I) TERMS IN ASCENDING ORDER. 

C 

K=NDAT A- 1 

IF(K.LE.O) GO TO 30 
DO 20 1 = 1, K 
N= NDAT A- I 
ISTOP= 0 
OO 10 J= 1 , N 

IF (X(J) ,LE.X(J+1)> GO TO 10 
SAVE=X ( J) 

X(J)=X(J+1) 

X( J«-1)=SAVF 
I3TOP=ISTOP+1 
1C CONTINUE 

IF(ISTOP.EQ.O) GO TO 30 
2 C CONTINUE 
C 

C SET CUMFRO(I) ARRAY 

C 

30 OO 40 1=1, NDAT A 
40 CUMFRQ(I)=I 
C 

C RESET SOME CUMFRO(I) ENTRYS IF X(I)=X(I+1) OCCURS 
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n 

00 50 1=2, NOATA 
50 CONI INvJt 

C PC ARE A = P(N) O p OBSERVATIONS 

00 759 1*1, NOATA 
759 PCAREA(I) = C'.JMFRQ ( I ) / D AT A 
PRINT 408, TITLE, SUBTITL 

408 F0RNAT<*1* ,55X ,2A1C,/,42X, ♦NORMAL DISTRIBUTION FITTEO TO ^,2 A10, 
1 ///) 

IF (RATIO, £0,0.0) GO TO 414 
PRINT 4 C2 , STRLEV, RATIO 

402 F0RMAT(29X , ♦STRESS LEVEL =*,F7.1,* »SI.»,16X 
1,14H STRESS RATIO =,F6.3//> 

GO TO 415 

414 PRINT 416, STRLEV 


416 F0RM4K29X, ♦STRESS LEVEL =*,F7.1,» PSI.*,16X 
1 , 23HST RESS RATIO = INFINITY//) 

415 PRINT 404, SUBTITL 
4Qh FORMAT (56X ,2A10,/» 

PRINT 403, (X(I) ,1=1, NOATA) 

403 FORMAT ( 6F2 1, 4) 

PRINT 3 

3 FORMAT (1H0) 

CALL MEAN 
CALL CHI SO A 
CALL DTE ST 
CALL ALPHA 

READ 502, IT, XLENGTH,Y LENGTH, X Ml, Y MIN, XMA ,Y MAX ,H LETTER 

502 FORMAT (Al, 9*, 7F10.0) 

READ 503, ROOT 

503 FORMAT (EA 10) 

IP (IT.NE. 1HN> GO TO 57 
CALL GRAPH 
53 AKURCY=. 30 001 
DO 54 1=1, NOATA 

NX ( I ) =(ALOG(X(I)/20.|+ALOG(20.) )MQQ00G. + .5 
X(‘l> = N X ( I) 

54 X(I) = X ( I) /10 0 003, 

PRINT 1, TITLE 

1 FORMAT (♦1*,55X,2A10) 

PRINT 401, SUBTITL 

401 FORMAT ( 33 X ,*LOG-NO.RMAL DISTRIBUTION FITTEO TO ♦,2AlO,///> 
IP (RATIO. EO. 0.0) GC TO 417 
PRINT 402, STRLEV, RATIO 
GO TO 418 

417 PRINT 416, STRLEV 
418 PRINT 2, SUBTITL 

2 PORMAT (49X ,»LOGS OP THE ♦,4A10,/) 

PRINT 413, (X(I) ,1 = 1, NOATA) 
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413 F0RMAT<6(3X,F12.4> ) 

PRINT 3 
_CA LL MEAN 

CAuTCHISaA' ' " 

CALL OTEST 
CALL ALPHA 

READ 5 0 2, IT »XLENGTH»YLENGTH»XMI,YMIN»XHA»Y‘iAX 
IFUT.NE.1HL> GO TO 57 
CALL GRAPH 
GO TO 7ifl 

56 CALL PLOT (0. ,0. ,999) 

CALL EXIT 

57 CALL PLOT (0.,0.,999> 

STOP 1111 

END 

SUBROUTINE MEAN 

C SUBROUTINE TO CALCULATE THE MEAN AND STANDARD DEVIATION OF DATA. 

COMMON CUHFROdOO) , NO AT A , X ( 1 0 0 > , OE V , XME A N, CLO ( 9 ) ,CU8<9) , FRF Q (9) ,K, 

1 CHI SQR, TITLE C2) , SURT ITL < 2) , CSV < 9 > ,CEV <9 ) , PC ARE A < 10 0 > , OSTAT(IQO) , 

2 AREA ( 9) , RE 0 A RE A (9) , EXFREQ (9) , U ( 9) , Z ( 10 0 ) * NX ( 1 j 0 ) * DA T A , AKURC Y t 

3 xm ax, xmin,psi,co,o,r,xlength,y length 
sigma= C.o 

DO 3 1*1, NOATA 
3 SIGMA=SIGMA+ X(I> 

XMEAN = SIGMA/OATA 
TOP2 = 0.0 
DO 9 1=1 , NOATA 

9 TOP 2 = TOP 2 ♦ C X C I ) - XMEAN)”2 

DEV =SQPT <TO?2/<9ATA - 1.0>> 

PRINT lu, XMEAN 


PRINT 15, DEV 

14 FORMAT (10X, ’SAMPLE MEAN =*,F14,4) 

15 FORMAT ( 10X ,’STD. OEV I AT I0N = * ♦ F12 . <♦) 

RETURN 
END 

FUNCTION PROB(X) 

THIS SUBROUTINE GIVES AREA UNDER NORMAL CURVE PROM -Z TO + Z 

WITH AN ACCURACY OF 0.00005 

Z VALUE GIVEN BY CALLING PPOGRAM MUST BE A POSITIVE NUMBER. 

IF (X-l. 2) 11, 11,12 

11 XSQ=X*X 

PRO 9= 0. 7978 S455*X*( 0.9999977 4- XSQ’l 0.1 665 943 3-XSQ* ( 0. 0 2 46 33 31 Q-XS 
10’C. 00239748fa7)> ) 

RETURN 

12 IF ( X-2 , 9 ) 13,14,14 

13 XSQ=X*X 
PR03=1. 0 
PT£RH=1. 0 
FACT OR= 1.0 
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ODDINT-3* 0 ' 

970 PTERM=-PT£RM*XSQ/<2.0*FACTOR) 

T£PM=PTERM/ODDINT 
PR09=PR0B+ TERM 

IF( A9S (TFRM) - 0.00007 ) 80,90,90 

90 FACTOR =FACTOR +1.0 

ODOINT = ODDINT4-2.0 
GO TO 970 

80 PRO9=0.79788495*X*PROB 

RETURN 

14 REC XSQ- 1.0 ✓ (X*X> 

PRO 9= 1.0 - 0.79788453*EXP(-X*X/2.0)/X* ( 1 . 0-REC XSQ* < 1. -R£C XSG* ( 3. 
1 - RECXSOM15. - RECXSQ*105. >>>) 

RETURN 

END 

SUBROUTINE CHISOA 

C SUBROUTINE TO FIT A HISTOGRAM TO THE DATA AND PERFORM THE CHI-SQUA 

C TEST FOR THE NORMAL OR LOG-NORMAL DISTRIBUTIONS. 

COMMON CUMFRQ(IOO) , NOAT A , X { 1 0 0 > , DEV , XME AN , CLB (9 ) ,CUB(9) ,FREQ(9> ,K, 

1 C HI SQR, T I TLE ( 2 ) ,SUBTITL(2) ,CSV( 9) ,CEV (9) , PC A RE A ( 10 0 ) ,DSTAT ( 10 0 ) , 

2 AREA (9) , REQAREA (9) , EXFREf) { 9 ) ,U(9) ,2(100) , NX (10 01 , Q A T A , AKURCY, 

3 XMAX.XMIN ,PSI,C0,0,R,XLENGTH,YLENGTH,YMAX,YMIN,XMA, XMI,H,COM(3) ,W 

! 01 MENSICN XFREO (9) 

i CHISQR= .0 

! n TO DETERMINE THE NUMBER OF CLASS INTERVALS, K 

K= 1 +3.3 *AL0G1Q COAT A) 

REALK=K 

;■ C IN OROEP’ TO DETERMINE THE RANGE, FIND X(MAX) AND X(MIN) 

XM A X = X ( 1 ) 

XM IN = X(l) 

DO 17 1=1, NOAT A 

IF ( xm.GT.XMAX ) XMAX = X(I) 

17 IF(X(I).LT. XM IN) XMIN=X(I) 

RANGE= XMAX- XMIN 

i c TO DETERMINE THE CLASS INTERVAL WIDTH, W 

: C ROUTINE TO ROUND OFF CLASS WIDTH TO SAME NUMBER OF PLACES AS THE A 

DIVIDE = l.C/AXURCY 

20 KW=( ( (RANGE+AKURCY) /REALK)+.5*AKURCY)*QlVlQt 
: R< 1 KW 

W = RK1/OIVIDE 
DO 22 1=1, K 
i A = I 

i B = 0 • 5* AKURCY 

i CSV(I)= XMIN+ (A-l. 0)*W 

j C£ V< I) = CSV(I) +W-AKURCY 

: C L B ( I ) = CS V ( I ) - B 

22 CU 9 ( I ) = CEV ( I ) + 9 
CEVCO = XMAX 
CU 9 ( K) = CEV(K) ♦B 
OO 23 J= 1 , K 
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XFREQ( J)=ti'.0 

23 FREQ(J)=0.0 

00 24 1 = 1, NDAT A 
00 24 J=1,K 

IP(X(I).G£.CL8CJ>.AN0.X(I) . IT.CUB(J) ) FREQ ( J) = FREQ ( J) + 1. 0 

24 CONTINUE 

no 25 J= 1 ,K 

25 XFR£Q( J)=FREO( J) 

Y=Q. 

PRINT 62 , X MAX 
PRINT 63 , XMIN 
PRINT 65, W 

C CHI-SQUARE TEST 

26 PRINT 41 
PRINT 4C6 
DO 3 0 1=1, K 

Z ( I ) = ( CUB ( I) - XMEAN) / DEV 
T= A 33 < Z<I> » 

3 0 ARE A ( I > = PROB (T) /2 • 0 

REQARZ A ( 1 ) = 0.5 - AREA (1) 

MANU=K-1 
DO 32 1=2, MANU 
M=I-1 

IF ( (Z(I) ,GE ,0,0, AND, Z (M) . GE . 0 . 0 ) . OR, ( Z < I) . tc . 0 . 0 . AN 0 . Z < M > ,l£ . 0 

1 ) ) GO TO 31 

REQAREA ( I > = AREA (I) +AREA <M) 

GO TO 32 

31 REOAREA ( I ) = ABS< AREAt I) -AREA <M> ) 

32 CONTINUE 

RcQAREA (K> = 0 . 5-AREA (K-l) 

00 90 M = 1 , K 

30 EXFREQ<M> =OATA*R£QAR£A(M> 

1 = 1 

2420 IF(FREQ(I) .GE.5.) GO TO 2430 

EXFREQ ( I + 1 ) -EXFREO ( I + l ) +EX FREQ ( I ) 

FREQ(I+l)=FREO(If 1) +FP£0<I) 

J = I 

00 2425 L = 1 , J 
EXFREO(L) =0. 

2425 FREQ (L ) = 0 . 

1 = 1 + 1 

GO TO 2420 
2430 1= K 

2440 IF (F RE 0(11 « GE . 5 . ) GO TO 2450 

EXFREQ ( I- 1 ) =EXF9£Q ( 1-1) +EXFREQ ( I ) 

FREQ(I-l)=FREO(I-l)4FPEO(I) 

00 2445 L= I, X 
EXFREQ (L>=0. 

2445 FRFQ<l)=0. 

1=1- 1_ 
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GO TO 244Q 
245 C CONTINUE 

00 85 M= 1 , K 
U(M) -0. C 

I c (EXF3£Q(m) « £Q . 0 . ) GO TO 85 

U(M) =( (EXPREQ(M) -FRED (M) )**2) /EXFREO(H) 

35 CONTINUE 
00 90 M=1,K 

9 D CHISOR=CHISQ R»U ( M ) 

0 -TO PRINT THE TA9LT FOR CHI -SQUARE' TEST " 

DO 33 1=1, K 

33 PRINT 34,I,CLB<I) ,CUB(I) , E XFREQ ( I) , FREQ ( I) » U( I > 

PRINT 35, CHISOR 
00 150 1=1, K 
150 FR£0(I)=XFRE0(I) 

62 FORMAT ( 1 OX ,*MAXIMUM V ALUE = * , F13. 4) 

63 FORMAT <10X*MINIMUM VAIUE=* , F13. 4) 

65 FORMAT (10X,*CLASS WIDTH=*, F15. 4) 

41 FORMAT (1HQ) 

406 FORMAT (5X,5H CE LL . 1 0 X , 1 OH LOWE R CELL,11X , 10HUPPER CELL , 1 3X , 3HEXPEC 

IT ED , 13 X ,8H OBSERVED, 13X, 1 1HCHI -SQUARE O/fiX , 6HNUMB£R,13X, 8H BO UN CARY, 
213X, 3HB0UN0ARY,13X,9HFRE0UENCY,12X,9HFRtQUENCY, 12X, 13HVALUE CF Ct 
3LL/) 

3* FORMAT (10X,T2, P 20,4,4F21.4) 

35 FOPMAT(*0*,77X,*TQTAL CHI-SQUAF.ED VALUE =*,Fl3.4> 

RETURN 

END 

SUBROUTINE DTEST 

C SUBROUTINE TO CALCULATE THE KOLMCGOROV-SMI RNOV Q-VALUIS. 

COMMON CUM FRO (10 0) , ND AT A , X ( 1 0 0 ) , OE V, XMt A N , CLO ( 9 ) ,CUB(9> , FREO (9) ,K, 

1 CHISQR, TITLE ( 2) , SUOTITL (2 > , CS V ( 9 ) ,CEV ( 9) , PC ARE A ( 10 0 ) , DS TA T ( 1 0 0 ) , 

2 A REA ( 9 ) , RE 0 A REA (9 ) ,EX FREQ (9) ,U<9>, 7(100), NX (100) , OAT A , AKU RCY , 

3 XHAX,XHIN,PSI,C0,0,R,XLEN6TH,YLENGTH,YMAX ,YMIN,XMA,XM1, 

4 M,COM(3) , W, ALPHA 3, ALPHA 4, FOOT (30) ,IT,SKS 
INTEGER TITLE, SU n TITL 

00 706 1=1, NDATA 

7(1) = (X(I) - XME AN) /OE V 

IF (7(D) 7 03, 704, 705 

703 T = A 3S ( Z ( I ) ) 

C ARNCN= AREA UNDER THE NORMAL CURVE TO LEFT OF Z FOR NEGATIVE 7. 

ARUNCN = (l.O-RPOB(T) J/2.0 
OSTAT(I) = A=>UNCN - PCAREA(I) 

GO TO 7D6 

704 DSTAT(I) = ,5 - PCAREA(I) 

GO TO 70 6 

705 T = 7.(1) 

C ARUNCP=A REA UNDER THE NORMAL CURVE TO LEFT OF 7 FOR POSITIVE 7, 

ARUNCP = DR09(T)/2.0 + .500 

DSTAT(I) = (APUNCP - PC AREA ( I > ) 

706 CONTINUF - 
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DO 90 I - 1 ♦ NDA T A 
90 OSTAT (I)=ABS(DSTAT (I) ) 

SKS-DST AT ( 1 ) 

00 100 1 = 2, NOATA 

IF ( DST A T (I ) . GT « S KS ) 3KS=DSTAT(I> 

100 CONTINUE / 

PRINT 70B.SUBTITL 

PRINT 707* (OSTAT (I) ,1=1, NOATA) 

707 FORMAT ( FF2C. 5) 

7 0B FORMAT <//40X,53H D VALUES FOR KOLMOGOROV-SMIRNOV GOODNESS QF FIT 
1TEST , / , 4i.x ,*LISTED IN THE SAME ORDER AS *,2A10,/) 

RETURN 

END 

SUBROUTINE ALPHA 

COMMON CUIFRa(lOC) ,NDATA,X < 10 0 > , OE V , XME AN , CL B < 9 ) ,CUB ( 9) , FRE 0 ( 9) ,K, 

1 CHISQR, TITLE ( 2) ,SUBTITL(2) ,CSV (9 ) ,CEV <9 ) , PC ARE A ( 13 0 ) , OSTAT ( 1C C ) , 

2 AREA (9) » RED A RE A (9) , EX FREQ (9) ,U ( 9) , Z < 1001 ,NX <10 0) .DATA, AKURCY, 

3 X MAX , X MIN, PSI,CD,D, RATIO, XLENGTH, YLENGTH ,YMAX, YMIN.XMA ,XMI , 

4 H,COM(3) ,W,ALPHA3,ALPHA4,FOOT(30) ,IT,SKS 

0 SUBROUTINE TO CALCULATE THE COEFFICIENTS OF SKEWNESS AND KURTOSIS 

C CALCULATE THE THIRD MOMENT OF THE DATA (SKEWNESS) 

T0 P 3 = 0.0 
VAR = 0.0 

00 710 I =1, NOATA 

VAR = VAR ♦ (XCI) - XME AN) **2 

710 TOP3 = T 0°3 + (X(I) - XMEAN>»*3 
SKEW = TOP 3 / DATA 

STOEV = SORT ( VAR/ DATA ) 

C AL PH A3 = MOMENT COEFFICIENT OF SKEWNESS. 

ALPHA3 = SKEW/ (STDEV**3) 

C CALCULATE THE FOURTH MOMENT OF THE DATA (KURTOSIS). 

TOP 4 = 0.0 

00 711 1-1, NOATA 

711 TOP4 = TO p 4 * (X ( I ) - XMEAN)**4 
TKURT = TOP4 / DATA 

C ALPHA4 = MOMENT COEFFICIENT OF KURTOSIS. 

AL PH A 4 = TKURT/ (STD EV**4) 

PRINT 712 
PRINT 713 

PRINT 714, ALPHA3, ALPHA4 

712 FORMAT ( / / /I 9X , 39HM0MENT COEFFICIENT OF SKEWNESS (ALPHA3) , 1 5X, 39HM 
10MENT COEFFICIENT OF KURTOSIS ( ALPHA 4) /) 

713 FO p MAT (21X,*FOR NORMAL DISTRIBUTION ALPHA 3 = 0. 0* ,21X , ♦FOR NORMAL 
1DISTR I BUTTON ALPHA*. = 3.0*,/) 

71*. FOPMAT ( 2BX , 25HFOR ABOVE DATA ALPHA 3 = Ffc . 3 , 2n>X , 2 5HFO R ABOVE DATA 

1 A L P H A 4 -, r 6 • 3) 

RETURN 

END 

SUBROUTINE GRAPH 

_ COMMO N CUMFPQ(IGO) , NOAT A , X ( 10 0 ) , DEV. XMEAN . CL 3 ( 9 ) .CUB (9) . FRc. H9) .K. 


o on ooo ooo oo 
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1 CHI$QR,TITLE(2) , SUBTITL ( 2 ) ", CSV ( 9) ,CE\/(9) ,°CaREA(10Q> ,DSTAT (ICO) , 

2 ARE A ( 9 ) , REQAREA (9) ,EXFREQ(9) , U( 9 ) , Z (10 0 ) , NX ( 1 0 0 ) , OAT A , A KU RC v , 

3 XMAX,XMIN,PSI,CD,D,R,XLENGTH,YL£NGTH,YMAX,YMIN, XMA ,XMI , 

4 H,C0M(3) , W,ALPHA3, ALPHA4,F00T (30) , IT,SKS,UNIT 
DETERMINE DEFAULTS GR SPECIFIED PARAMETERS 

IF(XLENGTH.EQ. 0. > XLENGTH=6. 

IF(YL£NGTH.EQ. 0. ) Y LENGTH-5 • 

IF<YL£NGTH.GT.5. ) YL£NGTH=5. 

DO 1 1=1, K 

IF(YMAX.LT.FREQd) ) YMAX = FRE0 ( I) 

1 CONTINUE 
I=YMAX 

IF( (1/2*2) .NS. I) YMAX=YMAX+ 1. 

6 XM IN = XMI N - 0.10 * (XMAX - XMIN) 

8 I p (XMIN . LT • 0.0) XMIN = 0.0 
XDIF=XMAX-XMIN 
H = 0.15 

DETERMINE SCALING FACTORS 

XSCALE = (XMAX - XMIN) / XLENGTH 
YSCALE = YMAX / YLENGTH 
E = IFIX(ALOG10(XMAX - XMIN)) 

STP = 10.0 ** E 

LOCATE PLOTTER PEN 

CALL PL0T(XLENGTH + 2. ,0. ,-3) 

CALL PLOTC C. ,-il. ,-3) 

CALL PLOT (0. ,5. ,-3> 

CONSTRUCT Y-AXIS 

CALL PLOT (0., YLENGTH, 2) 

DIV=1. /YSCALE 

21 IF(DIV.G£. (2.*H) ) GO TO 22 

DIV=2.*DIV 

GO TO 21 

22 ST £ 3 =DIV 
YY= 0 . 

23 CALL PLOT (.05,YV,3) 

CALL PLOT (-.05, YY ,2) 

YN=YY*VSC4LE 

CALL NUMDER ( - . 3 , Y Y , H , YN , 0 . , - 1 ) 

YY=YY+STEP 

IF (YY.LE . (YLENGTH+,01)) GO TO 23 
YY= ( YLENGTH- 3. 75) /2. 

IP (YY.LT.0.) YY=0. 

CALL SYM90L ( - . 4 , YY , H , 30HFREQUE NC Y/C LASS INTERVAL W I DTH , 90 . , 3 G ) 
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CONSTRUCT, DRAW, ANO LABEL X-AXIS 

CALL PLOT (0. ,0. ,3) 

IP(IT.cO.lHL) GO TO 25 
XMIN=IFIX(XMIN/STP) 

XM AX=I FIX (XMAX/ST P+1.0) 

XSCALE - ( XMAX - XMIN ) / XLENGTH * STP 

25 XDIF=XMAX-XMTN 
DIV=10.*XOIF/XLENGTH 

IF (XMIN.tQ.O. ) GO TC 26 
CALL PLOT (.3,0. ,2) 

CALL °LOT (.35,0. ,3) 

CALL PLOT (XLENGTH, 0. ,2) 

CALL PLOT (.35, .05, 3) 

CALL PLOT (.25,-. 05,2) 

CALL PLOT (.3, -.05, 3) 

CALL PLOT (.4, .05, 2) 

GO TO 28 

26 CALL PLOT < XLENGTH , 0 ., 2 ) 

28 IF(0IV.LT.12.7) GO TO 30 

DIV=DIV/10. 

GO TO 28 
30 OIV=DIV/10. 

I c ( DIV » LT. 0.2) OIV=9IV*10. 

32 CALL NUMBER ( 0 . , - . 2 , H , 0 . , 0 . , 0) 

XX- 0 . 

DO 35 1=1,25 
XX= XX+l./OIV 
IF (XX. GT. XLENGTH) 40,33 

33 CALL PLOT (XX,. C5, 3) 

CALL PLOT (XX,-. 05,2) 

XN= XMIN+ I *XOIF/(OIV*X LENGTH) 

IF ( IT. EC. 1HL) GO TO 37 

CALL NUMBER ( XX- . 1 , - . 2, H , XN , C . , 0 ) 

35 CONTINUE 

40 DO 4 1 1=1,2 

IF (SUBTITL (I) .EO.IH ) GO TO 42 

41 CONTINUE 
1 = 2 

GO TO 43 

37 CALL NUMBER ( XX- . 2 , - . 2, H, XN , 0. , 2) 

GO TO 35 

42 1=1-1 

43 XX= (XLENGTH-I) /2. 

1=1*10 

CALL SYMBOL ( XX , - . 5 , . 15 , SUBTI TL ( 1 ) , 0 . , I ) 
IF(IT.EO.IHL) GO TO 46 
IF(STP.LT. 1.01) GO TO 48 . 

CALL WHERE (XX, YY) 

CALL SYMBOL (XX,-.5,H,5H X 10,0. ,5) 
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"CALX ifflLP£"<XX;YY> 

CALL NUMBER ( X X+H , Y Y+ . 5*H, . 5*H , E , 0 . , -1 ) 

CONSTRUCT ANO DRAW HISTOGRAM 

GO TO 48 

46 DO *7 1=1, K 

CLB (I) = (CLB(I)-XMIN) /XSCALE 

47 CUB(I)=(CUP(T)-XMIN)/XSCALE 
GO TO 52 

48 DO 50 1=1, X 

CLB ( I) = ( <CLR(I>/STP)-XMIN)/(XSCALE/STP) 

5C CUR (I) = ( (CUB(I)/STP)-XMIN)/ (XSCALE/STP) 

52 CALL PLOT (CL3(1) , 0, , 3) 

00 55 1=1, K 

Y= FREQ (D/YSCALE 

CALL PLOT (CLB ( I) ,Y,2) 

CALL PLOT (CUB ( I ) ,Y,2) 

54 CALL PLOT (CUB ( I ) , 0 . , 2) 

55 CONTINUE 

COMPUTE AND DRAW NORMAL CURVE ONE POINT AT A TIME 

60 ST£°=XDI P *STP/10O, 

IF(IT.EC.iHL) STEP=STEP/STP 
C=l./(nev*2. 50665) 

XX=XMIN 

IP(IT.-:0.1HN> XX=XMIN*STP 
CALL PLOT ( 0# , 0 • , 3) 

DO 100 1=1,100 

Y=C*EXP(-»5MXX-XMEAM)* + 2/DEV**2)/YSCAL£*N0ATA*W 
XU= (XX-XMIN*STP) /XSCALE 
IF (IT. EG. 1HL) XIJ= (XX-XMIN) /XSCALE 
IF ( XU, GT . 2 0 . ) GO TO 10G 
IF(XMIN.EO.O. » 80,70 
70 IF(XU.GE.O .4) 80,90 
80 CALL PLOT (XU,Y,2) 

GO TO ICO 

90 CALL PLOT ( XU , Y , 3 ) 

100 XX=XX+STE° 

OTHER alpha-numeric COMMENTARY 

13C CALL PLOT ( 0 « » - 1 • * - 3) 

CALL SYMBOL ( 0. ,0. ,H, 11HMEAN VALUE* , 0. , 11) 

CALL SYM3OL(0 . ,-2.*H,H, 19HSTAN0AP0 DEVIATION 1 , 0. * 19) 

CALL SYMBOL ( 0 . , -4 . *H , H , 24H KOLMOG ORO V-SM IRNO V TEST*, 0., 24) 
CALL WHERE(XX.YY) 

XX = X X + H 

CALL SYMBOL ( 0 . , -6 , *H , H , 17 HCH I-S OU«REO TEST*, C. , 17) 

CALL SYMBOL (Q.,-8.*H,H,9HSKEWNESS*,0« ,9) 
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CALL SYMSOLCQ. , - 1 0 .*H, H , 9HKURT0S IS 1 1 0 . ,9) 

IFUT.tO.lHN) CALL NUMBER ( XX , 0 . , H , X ME AN , C . , 1) 

IF < IT. EO. 1HL) CALL NUMBER ( XX , 0 . , H , XMEAN , 0. , 3 ) 

CALL WHERE (XIN.YY) 

XIN=XIN+2.*H 

IF (IT. £0# 1HNI CALL NUMBER (XX , -2,*H ,H ,DEV , C . , 1 ) 

IFdT.EO.lHL> CALL NUMBER ( X X , - 2. *H ,H ,OEV , G . , 3) 

CALL NUMBER(XX,-4.*H,H,SKS, 0. ,3) 

CALL NUMBER(XX,-6.*H,H,CHIS0R,C. , 3) 

CALL NUMBER(XX,-8.*H,H,ALPHA3,0. ,3) 

CALL NUMBS R < XX, - 10. *H,H, ALPHA 4,0. ,3) 

CALL SYMBOL ( XIN, 0 ., H , UNIT , 0 101 

CALL .. SYMBOL_ < X IN , -2 . *H , H , UNIT , 0 . , 1 0) 

CALL SYMBOL ( 0 • » -14 • *H, H , FOOT < 1 ) » 0 . , 50) 

CALL SYMBOL(0. ,-16.*H,H,F00T(&> ,0. ,50) 

CALL SYMBOLCO. ,-18. *H,H, FOOT ( 11) , 0. ,5G> 

CALL SY*BOL<0. ,-20. *H,H, FOOT life) , 0.,50> 

CALL SYMBOL(0.,-22.*H,H,FOOT<21>,0.,50) 

CALL SYMBOLtO. ,-24.*H,H,R00T(26> , 0. ,50) 

CALL PLOT (XLENGTH+2, ,C. ,3) 

IF ( IT . EO . 1HL > GO TO 150 
XX=(XLENGTH-3.75) /2. 

IF ( XX. LT . 0 . ) XX = 0. 

CALL SYM30L ( XX, 6. 25 , H , 3CHNORMAL DISTRIBUTION °ARAMET£RS,0. ,30) 
RETURN 

150 XX=(XLENGTH-4.25)/2. 

IF(XX.LT.3.) XX=0. 

CALL SYMBOL ( X X , 6 . 25 , H , 34HL0G NORMAL DISTRIBUTION PARA ME T E RS , 0 • , 
1 34) 

RETURN 

END 
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APPENDIX D 
PROGRAM WEI BULL 

This program calculates the three Weibull distribution parameters 
(B, n, and y) from cycles -to -failure data. It uses those parameters to 
calculate cycle life for 99% and 90% reliability with a 90% confidence 
interval. It also performs the Chi -squared and Kolmogorov-Smimov 
goodness-of-fit tests. 

The program input consists of: 

1. A header card - to identify the data block. 

2. A set of data cards (in increasing order of cycles to failure) 

3. Trailer card - to separate data blocks. 


The input format for the header card is: 

a. Columns 2-7 date code in alphameric format 

b. 8-9 blank 


c. 


d. 


e . 

f. 


10-40 run identification in alphameric format 
41-46 number of data (sample size) in fixed point 
format 

48-52 minimum life increment in fixed point format 
53-80 not used 


Example: 

Dec 71 WEIBULL SL = 100,000 SR = Infinity 35 100 

The input format for the data cards is: 

a. Columns 2-7 cycles to failure in fixed point format 

b. 9-15 component life in fixed point format (same as 

for individual fatigue life tests) . 
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c. 20-23 date code (prints out only) in alphameric format 

d. 25-27 number of suspensions in fixed point format (either 

failures just removed from test, or remaining good 
items when testing by groups) . 

e. 28-31 blank. 

32-35 test code in alphameric format (Note: 28-42 simply 
prints as punched and can all be left blank or used 
as comments space) . 

36-39 lot code in alphameric format. 

40-43 blank. 

44-80 not used. 

The input format for the trailer card is: 
a. Columns 1-8 not used. 

9-12 punched: - 1.0 in floating point format. 

13-80 not used. 

By use of the trailer card, the program is written to process as many 
sets of data as desired. An end of file card is used to signal the end of 
the last data deck. 

Following is a list of important variables and symbols used in 


program WEIBULL: 



List of Definitions for Program WEIBULL 
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Main Program: 

NDAT = number of data points 

W = NDATA + NODATA = DATA =W=V=J=K=FNZ= NODAT 
IACTL = cycles of test at which a failure occurs 
IX»Z = cycles of life of specimen at failure 

MINL * IACTL 

INCR = minimum life increment 

C * cumulative life increment 

R = median rank 

Rj * log (log) median rank 

Qj_ B log CX - G) 

* B = weibull slope, 3 
BETA = inverse of slope 

Ej * E = goodness of fit of regression line 
Gj * minimum life parameter, y 

U = median life 

ETA * scale parameter ■ n 

CG = plus confidence interval 

ONE = 1% life = 99% reliability 

CONE = upper confidence limit on 1% life 

CONM » lower confidence limit on 1% life 

TEN '= 10% life = 90% reliability 

CTEN = upper confidence limit on 10% life 

CTENM = lower confidence limit on 10% life 

CU = upper confidence limit on median life (U) 

CUM = lower confidence limit on median life 



Subroutine for Kolmogorov-Smimov Test 
(Subroutine DTEST) 
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Z(I) = 
CUMFREQ(I) = 
PERCF(I) 
DSTAT(I) 


Subroutine 

K 

XMAX 

XMIN 

CSV 

CEV 

CLB 

CUB 

FREQ(J) 

REQAREA (J) = 

EXFREQ(J) 


Weibull cumulative frequency distribution 
cumulative observations (numher of failures) 
data cumulative frequency as percentage of failures 
absolute difference between the data cumulative frequency 
and the hypothesized cumulative frequency 
for Chi -squared test (Subroutines CHISQA and CHISQB) 
number of cells 

largest value of cycles to failure 

smallest value of cycles to failure 

cell starting value 

cell end value 

cell lower bound 

cell upper bound 

number of observations in J** 1 cell 
expected value of J** 1 cell (percentage) 
expected number of observations 


The program WEIBULL listing in Fortran language follows : 
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C WEI BULL 

C—-—- WEI BULL 

C -WEI BULL 

PROGRAM WEI0ULL ( INPUT, OUTPUT, TAPE2aINPUT,TAPE3aOUTPUT, 

1 TAPEXilARE^lAREStPLOTsTAPEl) 

DIMENSION y<100) ,FREQ(9 ) ,Cl 9(9) »CUB (9) ,OSTAT(100) tCUMFRQ(iOO) , 
1 PERCF(IOO) 

CQMMQN/N4/ Is P . 

COMMON /N5/ IACTL* IX* Al* A2, NOR* A3* A4* A5* A6 

1 J ■ 0 

. K..1 Jl - 

I a 0 
ICOUN * 0 

NOq. a 0 

XPREV ® 0*0 
RPREV a 0*0 

REWIND A ... 

REWIND 5 

READ (2*10) NODAT* INCR»M 

IF_tEuF (2j ) 03 *j8q. 

I 10 FORMAT (40H ,216.18) 

80 IF (NODAT-50) 82*83.83 

02 WRITE (3 iM) 

40 FORMAT (5H1DATE.15X.9HRUN IDENT ) 

WRITE (3*10) 

; 7 READ (2 »201 IACTL-* IX J Al» NOR «■ A3* A4i AS* A6 

20 FORMAT (I* * IX* 17, lx, 2A4, l3, 4A4, F16.6) 

IF (IX) 4, 3, 2 

2 J a J ♦ 1- . 

| 3 K a K ♦ 1 

! I a X ♦ 1 

! Y < I X ■ IACTI _ ... - - • 

IF (NOR) 6, 6, 5 
6 NOR a 1 

5 NOD a. NOD ♦ NOR .... 

WRITE (4) IACTL, IX* Al, A2* NOR, A3, A4, A5, A6 
60 TO 7 

A IF (NODAT - NOD) B, 9* 8 . 

C— — ERROR STOP - NO, OF DATA NOT CORRECT 


I 8 WRITE (3*30) NODAT, NOD ! 

30 FORMAT (22H1N0. OF OATA INCORRECT /1HO* 16 , 5x» 16 / 1H1) ! 

. RFWTnD 4 ... : 

60 TO 83 j 

9 REWIND 4 | 

WRITE (3 .220) NODA T. J. IN CRlM j 

220 FORMAT (1HO»4x,UHNO. OF DATA I6,10X,13H NO. OF FAIL I6.10X*! 
1 15H MIN LIFE I NCR I6*5X,5H M s .16//) | 

un f j 1 a i 5 
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230 FORMAT ( 2X »6h ACTUAL » 3x » *HCOMP # 56 x » 4HC0MP/ 

13X.4rtLIFE ,4X,4HLIFE ,9H P DATE ,4H N0R.6X, 10HTEST LOT 
... _.21LX»11HMEC1AN..RAHK.4 8X*4HLIFE- //) 

Z - 0,0 

W a NODAT 

NDA T A a NOOAT ' “ 

■MOQAIA-g- NOD AT - . . ... 

DATA * NOOAT 
AKURCY a l, 

...V JL W 

READ (4) IACTL* IX. Al. a2» NOR. A3, A4, A5, A6 
REWIND 4 

| JHML..-S-.IACJL- - 

I FNZ a J 

j IF (MINL - 1) 11. 11, 12 

i .. U-._MlNL--a.-l - . - - 

I NCR a 1 

GO TO 13 

1 2 MlNL aM INL" l — _ 

| 13 Rl a 1,0 - 2,**(-l,0/W) ♦ (1.0 - 2. ##(1,0 - l . 0/W ) ) / ( W - 1.0) 

REWIND 4 

! DO. 14 — I _ a — 1*. K. 

! READ (4) IACTL, IX, Al, A2, NOR. A3, A4, A5, A6 

! X a IX 

IEU) 16 » 16* U- — - 

16 IF < I * 1) 18, 18, 19 

19 IF(XPREV) 18, 18 , 21 

21 &I.M- RPR£Y . 

18 K 1 a NOR 

00 22 «2a 1, Kl 

. V a y ■ 1 , Q 

RINV«l,0*(W-l,0)*(Rl-1.0*2,0**(-l,0/W) ) / (2,0** ( 1 . 0-1 . 0/W ) -l .0 ) 

22 Z a WlNV ♦ (W ♦ 1.0 - RINV)/(1,0 ♦ V) 

R a...O^. ... 

GO TO 15 

17 IF ( I - 1) 24, 24, 23 

23 .-jF-lAgflEVl £8 j....26 > .._24 

24 Z » i . (W ♦ 1,0 - Z)/(1,0 ♦ V) 

26 Ral,o-2.0**(-1.0/W)*( (Z-1,0)/ ( W-l .0 ) ) * ( 2. 0** ( 1 . 0-1 . 0/W) -1 .0 ) 

V- a ¥. _a_.ua . .. ..... .. 

15 XPR£v » X 
RPREV a R 

WR IT E(3 ,241 ).--1AC-TLj IX, A.U A£*-NQft. A3. A4.,. A5, A6.R.X 

241 FORMAT(1X,I7,1X.I7,1X,2A4,1X,I3,3X»4A4.F16,6,F12.0) 

ICOUn a ICOUn ♦ 1 

IF.. ..(ICQUN - 6EJ 14, Alt 41 ... 

41 ICOUN a 0 

WRITE (3,400) 

— A&O FORMAT- (1H1) . 

14 WRITE (5) X, R 
REWIND 4 
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! REWIND 5. ...... - 

El * 0.0 

IF ( Y ( 1 ) ,LE. 1000) GO TO 117 

IF { Y { 1 ) .ST.«. -UUUL.«ANO*- UU- *LE. LQjQlQQ) GO TO 116 
IF(Y(1) ,GT. 10000) INCR ■ 10000 
00 TO 127 

116. INCR ■ LOGO.- - - 

00 TO 127 

117 INCR a 100 

1 IF (M ,GT, 1000) M « M - 1000 

U7 El a Q«G-. - - ----- — - - 

DO 27 KN«*M»MINl*INCR 
| QSUM a 0*0 

.RSUM s.Q^.Q . . ..... — 

QSOS a 0,0 
RSOS = 0,0 

1 .... . PRQD.a.JM ,. . 

; G * KN - l 

: DO 26 I s 1, K 

READ ( 5) X tR - 

IF ( X) 28, 28, 29 
29 IF C X - 6) 31* 31* 32 

C-^ NEGATIVE ARGUMENT IN LOG. FUNCTION 
C— — 

. 21 jyRij£-...ta 4 lMi.. — 

100 FORMAT ( 1H1.26HNEG, ARG. IN LOG, FUNCTION /,1H1) 

GO TO 83 

. 32 - --O-l-a-ALGG-U -« G> 

IF (R ) 31, 31, 34 


mmmwm. .. NEGATIVE ARGUMENT IN LOG -FUNCTION 


34 Rla aLOG(ALOG(1,0/(1,0 - R) ) ) 

. . - -QSU M s QSUM » -Ql ... 

RSUM a R$UM ♦ Rl 
QSOS a QSOS ♦ Ql * Ql 

. RSOS s RSQS R.1.JR-4U 

PROD a PROD ♦ Ql * Rl 
28 CONTINUE 

B a (FN2 • PROD • QSUM * RSUM)/(FNZ * QSOS - QSUM • QSUM) 

A a (RSUM - B • QSUM) / FN2 

_C _.a_iiQRT 11EN2. f.JlSQS. - QSUM. * QSUM).. 9 (FN2 * BSOS-* RSUM * RSUM) 

E a <FNZ*PROD - QSUM * RSUM)/C 
IF CE - El) 37, 36, 36 

Gl a G 

Bl « B 
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' 27 CONTINUE 


37 

19. 

38 


M a G ♦ 1 • INCR 

IF ( i mCR - 1000.1 3a»UTti 16 

ARG a EXP ( -Al) 

BETA a 1.Q/B1 

ETA » LQ^i99A7*ABaJ **B£TA 

U a G1 ♦ (0*6931 5* ARG) ••BETA 
ONE a G1 ♦ (0.01005»ARG)**BETA 
I£N_s 


B a Hi * U.O ♦ 1.163 / SQRT(W)) . . , A _ 

CGaGl*.5*(U-Gl)«(4.32159**(1.0/B)-0.074**(l.0/B) )✓(«••( l .0/B) ) 
CO NST a Xl*641 JLJJi—/. JSQRT-(J4---f. 0*A91l5- . £*. BETA) 

FACT = CONST • 10.010030 * <0.995E-02*»8ETA/81 ) 

CONE a O NE ♦ FACT 


CONM a ONE - FACT 


CTEN a TEN ♦ FACT 
CTENM a TEN -'FACT 

EAC-t-JL CONST #Q.693l5**BETA/(Bl *0.691154 

CU a U ♦ FACT 
CUM a u - FACT 

WRITS— (3.450)- Bl. El < ETA-. - &1 .-CB. -BNE*~ CONE.-- -CONM 

WRITE (3»45n TEN. CTEN. CTENM. U. CU. CUM 
!450 FORMAT (// 12X, 13HWEIBULL SLOPE. 5X, 15H0OODNESS OF FIT ,5X. 

-11 5HSCALE. PARAMETER »/5X»3F2q. 5. //_ 

213X.12HMINIMUM LIFE .11X.9HPLUS CONF./ 

| 35X.2F20.S.// 

! - 4 . 9X a j teNQ N E- P ERCE NT LI F E sll4(j_9 H etUS--CON£- -.-14UU14MMLNUS CONF 

55X.3F20.5 //) 

451 FORMAT (1H*,9X,16HTEN PERCENT LIFE. 11X.9HPLUS CONF.IOX. 
.- H IM M INUS CONF -*y5JU 3F2 0.5.// ..... . . 


i 214X ♦ 1 1HMEOIAN LIFE » 1 IX .9HPLUS CONF *10X» 10HMINUS CONF ./ 

| 35X.3F20.5 //) 

L CAL) — DI ES ! (Y. B 1 »G 1 .U. NODATA^DS-T A T^PEaCF^C UMFRQ .E-TA. SKST AT ) 

j CALL CHISQa ( Y. DATA » ndata.prob.akurcy.x MEAN. DEV, Z.G l *B1 .ETa.FREQ 
1 XMAX.XMIN.CLB.CUB.NUMINTS.CELLWO) 

CA L L C UlSfi B (X.^lOATA.Gi,Bl>ETA^CMlSqR) 

CALL GRAPH (FReO.XMAX.XMIN,CLB.CUB.SkSTaT,chISQR,bi,eTa,G1.NOaTa 
1 NUMINTS.CELLWD) 

; GO TO 1 . 

! 83 CALL PLOT (0.0. 0.0. 999) 

; call exit 

L END . 



r> o r> ooo o oo 
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SUBROUTINE DTEST ( Y» B1 «G1 *U» NODATA » DS TAT »PERCF, CUMFRQ, ETA » SKSTAT ) 
DIMENSION Z(100) • Y (100) »DSTAT ( 100 ) , PERCF ( 1 00 ) , CUMFRQ ( 1 00 ) 

SUBROUTINE TO CALCULATE THE KOLMOGOROV-SMIRNOW O-VALUES 

DO 500 1 = 1 tNODATA 

Z(I> 8 1,0 - EXP (-( ( (Y(I)-Gl)/ETA)**Bl) ) 

500 CONTINUE 

SET CUMFRQ (2) ARRAY 

DO 501 I 8 1 tNODATA 

501 CUMFRQ (I) 8 I 

PERCF s>(N) OF THE NUMBER OF DATA 

DO 502 I si » NOD AT A 

502 PERCF (I) 8 CUMFRQ(I)/NODATA 
DO 503 Isl , NODAT A 

503 DSTAT 111-- 8 2(1) - PE*CF(I> 

PRINT 520 

PRINT 521 » (DSTAT (I) » 1 = 1 tNODATA) 

521 FORMAT lfe(lUX*FlO*5)> 

520 FORMAT (//40X.53H D VALUES FOR KOLMOGOROV-SMIRNOV GOODNESS OF FIT 
lTEST/41X»52H (LISTED IN THE SAME ORDER AS CYCLES-TO-FAILURE DATA)/) 
SKSTAT 8 --0«0 • 

DO 10 I»1 tNODATA 

IF(ABS(DSTAT(I) ) .GT.SKSTAT) SKSTAT 8 A8S(DSTAT(I) ) I 

10 CONTINUE... | 

PRINT 400 » SKSTAT I 

400 FORMAT(//»lOX»*KOLMOGOROV-SMlRNOV TEST RESULT 8 *,FR.5,/J j 

Return 

i 
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SUBROUTINE CHI SQA(X, DATA, NDATA, PROS, AKURCY.XMEAN, DEV, Z 
1 * Q.1 j-WJU ElAjf RE£L* XMAUaXMINj CLB • CuB ,J4, * ) 

c — — -subroutine to fit a histogram to the data and perform the chi-squa 

|c test FOR THE WEI8ULL DISTRIBUTION 

[ DIMENSION X(5Qlj. CSVm.t_ CEV{9) * CLB (9) « _CUB491* 

IREQArEA (9) #AREA(9) ,EXFREQ<9) , FREQ(9), U<9) 

CHlSURa .0 

jC-----TO DETERMINE THE NUMBER . OF . CLASS INTERVALS,*- - 
j K a 1.0 ♦ 3,3#ALOG10 (DATA) 

j REALKsK 

jC- -IN ORDER. TQ DETERMINE. THE RANGE.*FIND X (MAX.) ANO X4MXN-I . 

! XMAX«X(1) 

{ XMlNs X(!) 

1 DO 17. I»1 fNDATA 

j IF ( X(I>.OT.XMAX ) XMAX a X ( I > 

;17 IF (X (I )/*LT • XMJN) XMINaX(I) 

| RANGES X.MAA- XMIN _ 

$C TO DETERMINE THE CLASS INTERVAL WIDTH,* 

iC ROUTINE TO ROUND OFF CLASS WIDTH TO SAME NUMBER OF PLACES AS THE A 

DIVIOE 8 I « 0/AKURC Y . 

KW 8 ( { (RANGE*AKURCT)/REALK) ♦*5#AKURCY) *oivide 
RK 1 a KW 

W « RK I. /DIVIDE . 

! PRINT 141 

PRINT 241 

| PRINT -_I77-*-NQATA, GX» Bl« ETA 

PRINT 4 1 
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IF (I ,EQ. 1) 00 TO 51 

IF (I .GT. 1 .AND, I ,LT. K ) 60 TO 52 

REQABEAiKi— ■ 1.0 AREA (K) 

GO TO 30 

REQAREA ( I ) a AREA (I) • 


GO TO 30 

-RE^aKEAXI) — a..AR£A.< I )_ ■ A R EAt-W) 

continue 

DO 80 M s 1,K 

EXFR EQ (M l «DATA» REQAg£A (M) 

U (M) ■ ( ( EXFRE0{M)«FREQ(M))##2)/EXFRe0(M) 

CHISuR*CHISQR*U(M) 

•-T.Q giOJV.T....THE. TABLE-FOR.-CJAI«SQUAR£-T£SJ- - 

DO 33 I s 1,K 

PRINT 34,ItCuB(I) *CU8(I)* EXFREQ(I) * FREQ (I) *U(I) 

-PRINT- 35-,~CHISQR - 

FORMAT ( lOXt 14HMAXIMUM VALUE«*F15.6) 

FORMAT < 1QX* 14HMINIMUM VALUE** Fl5,6) 

-XORMaTL-IOX, 12HCLASS BIDTh** F-17-.A) 

FORMAT (1H0) 

TFD?i 3X !f HORSFRUED » 1 | 1 SrSt ^qflABPft * a* * * °HUPPER CELL*13 x,8HEXPEC 

213X. 8HB0UNDARY*13X*9HFRE0UENCY*12X,9HFRE0UENCY»12X,13HVALUE Of CE 
3LL/) 


FORMAT (1HO*81X.25HTOTAL CHI-SQUARED VALUE s 
FORMAT (1H0»70X,#CHI«SQUARED TEST#,//) 

FO R M AT U H Q i 7 QX* »FIX£D CEU. - B 1 07MSA*// ) 

FORMAT (1H0, 10X* #INPUTS « **110* 3F15.3*/) 
CONTINUE 

REJ.UR.N 

END 


iFIO.B//) 
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SUBROUTINE CHISQB ( X * ND AT A , G1 » B 1 » ET A * CH I SOR ) 

QIMEtiS-LQN X(50)t CLB(9J» CUB (.9J-»_ REQAREA-UJ > AREA(9>» FREQ (9 ) « 

2EXFREQ (9) i U(9) 

—subroutine to fit a histogram to the DATA And perform the 
— CMI«SQUAB£D-T£ST .FOB THE IB ULL -D-LSTRIBUTION 

PRINT 341 


CHI SOR a 0,0 

. J a 0 . 

00 2b K*5»NDATAt5 

J 9 J* 1 

_IE . IK-U.T-* N0ATA1 CUB ( J) X (K) *-X (X-*-i)iy2-,Q ... 

IF (K ,EQ. 5) CLB(J)»X(1) 

IF (K ,GT» 5 , AND * K ,LT • NOATA) CUB (J)aCUB(J-l) 
— L-a_ 4 N 0 AIA«jq - — . 


IF(U.NE.O) AREA (J) al, O-EXP (-( ( (CUB ( J) -Ql ) /ETA) «M>B1 ) ) 
IF(L.EQ.O) AREA(J) >1,0 

FREQ ( U) 0 — - 

IF(J.EQ.l) REQAREA ( J) aAREA ( J) 

IF(J.GT.l.AND.L.NE.O) REQAREA (J> -AREA (J) -AREA (J-l) 


GO TO 26 

27 IF(U.NE.O) J«J*1 

CUB(j)aCuB(J-l) 

REQAREA (J) *1. O-AREA (J-l) 

IF < -( .- -^£♦-04 FR&Q 4U4 -» U 

26 CONTINUE 


I ■ J 

no 2 S J al«I . 

EXFREQ(J) a N0ATA*REQAREA(J) 
25 CONTINUE 


K-a~4-- ... 

62 I a 1 

; 2420 IF (EXFREQ(I) ,GE, 5.) GO TO 2430 

- .EX£R£QU*-1) » -EXFRECttI ♦ 1 ) * EXFREQ (I ) 
FREQ ( I ♦ I ) a FREQ (1*1) ♦ FREQ(I) 

J a I 

l QQ_ 24 .25._L5l.».J . . .. 

j EXFREQ(U) a 0, 

! 2425 FREQ(U) a 0, 

; . . 

GO TO 2420 
2430 I a k 

2440 ip (EXFREQ (I ) ,G£, 5,) GO TO 2450 

EXFREQ ( I»1 ) a EXFREQ ( 1-1 ) ♦EXFREQ ( I > 
FREQ(I-l) * FREQ(I-l) ♦ FREQ ( I ) 

DO 2445 4.^1 ,K 

EXFREQ(L) ■ 0, 

2445 FREQ (L) = 0, 
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L I -a. XjbI 

' GO TO 2440 

I 2450 CONTINUE 

!- - D0-8S Mal»K- . . 

U (M) s 0.0 

IF (EXFREQ(M) ,EQ, 0.) GO TO 85 

U(M) a ( (EXFREQ (M) - FREQ (M) ) «#2/EXFREQ (M) ) 


CHISQR a 0,0 
DO 90 Mai,K 

90 — CHISQR ■ CHISQR * . U(M) - 

J • K 

88 DO 3J 1 s hJ 

.33 PRINT 34 > UCL&(I-l»CUB (I)» EXER£Q( I) ,.ER£Q-< I-ULU-I > 

PRINT 35* CHISQR 

34 FORMAT (10X*I2.5F2l,6) 

35 F - O R m a T — U ilO»8aX»25RTQlAL- CH I - SQUARES V-AtUE--»*F10 .6) 

341 FORMaT(1H0»65X**VARIABLE CELL WIDTHS**//) 

78 CONTINUE 

... RETURN 

END 


ooo ooo ooo 


173 


SUBROUTINE. QRARH (FREGiXMAX «.XMI N i CLB* CUfl-*SKS • CH I SQR « BET A • ET A * 

1 GAMMA, NDATA,K,CELLWD) 

DIMENSION FREQ (9) *CLB (9) »CuB (9) 

INTEGER EOQTI3.0) •SUBJITL-U) 

LOGICAL NITIAL 

DATA NITIAL / .TRUE. / 

READ PLOT CARD 

.READ-l24.40Xn IPLOT^SUBTIIi.-* I PEN 
AOO FORMAT (4A10) 

IF (EOF (2) ) 999,3 

3 IF ( IPLQT»NE»7.HREIBUEL) GQ TO 999 
IF(NITIAL) *,5 

4 NITIAL ■ .FALSE. 

PEN. a 0.3 

IF ( IPEN .EQ. 10HBALL POINT) PEN a 0,0 
CALL INITIAL (0,I»PEN,0) 

5 READ 12*401 ) FOOT. . 

401 FORMAT (5Al0 ) 

DETERMINE DISTANCES 

2 XLENGTH * 6»0 
YLENUlM 0 

H a 0,15 
YMAX a 0,0 

UO l_Ial»K 

IF (YmAX.LT .FREQ ( I ) ) YMAXsFrEQ ( I ) 

1 CONTINUE 

FRINGE <UA0*NQATA -- 
FRINL a 0.25*NDATA 

IF(YMAX .GE. FRINL .AND. YMAX ,LT. FRINGE) GO TO 7 
IF (YMAX .GE. FRINGE) YMAX a IFIX { YMAX ♦ FRINGE) 

GO TO 6 

7 YMAX * IFIX(YMAX ♦ FRINL) 

6. CONTINUE _. . 

I a YMAX 

IF( (1/2*2) .NE.I) YMAXaYMAX*!, 

XMIN a X.MIN_- 0.10 * (XMAX - XMIN) 

XOIF a XMAX - XMIN 

DETERMINE. SCALING FACTORS 

XSCALE a (XMAX - XMIN) / XLENGTH 
... . YSCALE. «- YMAX-./...YLEMGTH .. 

E a IFIX(ALOGlO (XMAX - XMIN)) 

STP a 10,0 ** E 


o o o 


locate plotter pen 

CALL PLOT LXLENGTH # 0 • t -3 ) 
CALL PLOT (9 • **1 1 • t**3) 

£Akk £ kOI ( O ft 5 i ) 
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C 

jC 


CONSTRUCT Y-AXIS 


£ 

1C 

\C 


CALL PLOT ( 0 • * YLENGTH*2) 

DIV »1 . / Y SCALE 

21 IF(DIV,GE. (2.*H) > GO TO 22 
DIV*2.*DIV 

GO TO 21- - - 

22 STEP»DIV 
YYaO. 

23 CALL -PLOT (.Q5,YY,3) 

CALL PLOT (-.05,YY,2) 

YNaYY*YSCALE 

CALL-NUMBLS. . ( « . 3 « Y Y * H^ YbU 0 •-* »13 - . 

YY»YY«-STEP 

IF(YY.LE. (YLENGThU.Ol) ) GO TO 23 
Y Y.*J YLENGTH-3* 75 ) /2. 

IF (YY.LT.O.) YY«0, 

CALL SYMBOL (-o*, YY,H,30HFREQUENCY/CLASS INTERVAL WIDTH, 90, 


3 0)! 


CONSTRUCT t DRAW, AND LABEL X-AXIS 

— - CA LL P LOT- -UX» ,0 ,43) 

XMlNnIFIX(XMlN/STP) 

XMAX a IFIX (XMAX / STP * 1,0) 

- -XOIF a XMAX - XMIN . 

XSCALE a XDIF / XLENGTH * STP 

div«io.*xdif/xlength 
IFXXJI tN^EQ^D-* ) GO TO 26 
CALL PLOT (.3*0. *2) 

CALL PLOT (.35,0. ,3) 

- CALL-BLOT -(XLENQTH»0. »3) - 
CALL PLOT (.35, .05,3) 

CALL PLOT (.25, -.05, 2) 

. CALL- BLQ-T -U3* -.05*3) 

CALL PLOT (.A, .05,2) 

GO TO 28 

26 C ALL p L fi T (X L E NGTH, 0 . ,2) 

28 IF (Dl V »LT • l2.T ) GO TO 30 
DIV.OIV/IO, 

GO. TO 2fi 
DIV»OIV/10, 

IF (OIV.LI.0.2) DIV«DIV*10. 

2, H , 0 .,0. vO ) 

XX«0. 

DO 35 Ialt25 


30 


IF (XX. GT. XLENGTH) 40,33 
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33 CALL PLOT (XX».05»3) 

CALL PLOT (XX.-.05.2) 

XNaXMlN*I*XDIF/(DIV*XLENGTH) 

CALL NUMBER (XX-.1*-.2*H*XN*0.*0> 

35 CONTINUE 
40 00 41 1*1*2 

IF ( SUB TITL(I) ,EQ,1H ) GO TO *2 

41 CONTINUE 
U2 

GO TO A3 

42 1*1-1 

43 XXaTJtLENGTH»I)/2, 

1*1*10 

CALL SYMBOL (XX*-.5 # .l5,SUBTITL(l) *0.*I) 

IF (SIP..LI .-1 • Oi ) GO TO AS 
CALL WHERE (XX. YY. IFAKE) 

CALL SYMBOL (XX*-.5*H*5H X 10*0. *5) 

CALL-WHERE XXX* YY» IFAKE) 

CALL NUMBER (XX*H.YY«*,5*h.,5*h.E.0..-1) 

C 

. C COMST.RUCI~.AND. dram, hisxo.gr am 

c 

48 DO 50 1*1 *K 

CLBXii *XXCLBXl)VSTPj -XMIN) / (XSCALE/STP-) 

50 CUB(I)»( (CUB(I)/STP)-XMIN)/(XSCALE/STP) 

52 CALL PLOT (CLB ( 1 ) .0. *3) 

DO 55 1 *1 *K 

YbFREQ(I)/YSCALE 
CALL PLOT (CLB(I)*Y*2) 

CALL PLQT-XCU&U) .Y*2)-. _ 

54 CALL PLOT (CUB ( I ) , 0 . . 2) 

55 CONTINUE 

C .... 

C COMPUTE AND DRAW NORMAL CURVE ONE POINT AT A TIME 
C 

60 STEP._* XDIF *.STP / .150.0 
XX * XMIN * STP 
CALL PLOT (0 • • 0 • * 3) 

EACT -B-£XP_(-BETA) __ 

CHK ■ EXP(-l.O) 

IF (BETA .LT. 1.0) YSCALE*YSCAL£*EACT/CHK 

BE «-8£lA~Y_ ETA. ... . 

IDOIT ■ 2HNO 
DO 100 I s 1.150 

IF (GAMMA .GT . XX) QO T O 100 . 

Z » (XX - GAMMA) / ETA 

_Y * BE * Z •• (BETA - 1.0) * EXP ( - 1.0 • Z **_ BETA) 
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_ jL Y,SCAii:;OMTA. * CELLJrfQ . 

XU« (XX-XMIN«STP) /XSCALE 
IFdUOIT .EQ. 2HN0) GO TO 90 

IF(XU.G T.2Q. ) QO TQ 1 00 

IF(XMIN.EQ,0.) 80*70 
y 0 IF(XU.GE.O.A) 8o»’0 

.8.0 CALL. PLOT. (XUi.Y 1.2) 

GO TO 100 

90 caul plot <xu*y, 3) 

IDQLT»._3HYES | 

100 XX»XX4STEP 

C i 

C. _QIHE8... ALPHA-NUMERIC COMMENTARY j 

C ' ' 

130 CALL PLOT < 0 . * • 1 . * -3 > j 

CALL SYMBOL (0 . 0*0 *6 *H * 2 4 H K 0 L M 6 G 6 R 0 V - 5 MIR NO V TEST t To .6 *24 ) ' ~j 

CALL- WHERE- jAX.» YY 0 l FAKE) — - 

XXsXX+H 

, CALL SYMBOL (0.0,-2.0#H,H,l7HCHI-SQUARED TEST*. 0.0, 17) 

CALL-_S.YM.BQL-. 10-4 0 » "4.* Q M *2 1 HWEXillLL-SLOPE ( BETA-) 1.0.0.21) i 

CALL SYMBOL (0.0.-6.0*H,H,21HMINIMUM LIFE (GAMMA ) I ♦ 6 . 0 . 21 ) 

CALL SYMBOL ( 0 • 0 * -B. 0»M , M,a2MSCALE PARAMETER < ETA ) * 1 0 . 0 , 22) 
CALL-.miMBEa-UU*0.04H^SKS*11.0->3)- - . 

| CALL NUMBER (XX.-2.0*H,H,CHISQR,0.0,3) 

j CALL NUMBER (XX*-4.0*H»H*8ETA*0.0*3) 

CALL -NUMBER. (XX* *6.0#H*H+GAMMA^0 .0-*.«l.) . 

CALL NUMBER (XX.-8,0«H,H 9 EtA90.0*-1) 

CALL SYMBOL (0., -14. *H,H, FOOT (1 ) ,0., 50) 

! CALL S YM B O L (-0-.-*«l^.-*H *H , F OQT|^*L^SOL 

! CALL SYMBOL(0..-l8.*H,H,FOOT(ii) *0*.50) 

! CALL SYMBOL(0..-20.*M,H,FOOT(16) *0.*50) 

I CALL^SYMBOLIO.. . -E2. »M^M4FC0T ( 2 1 ) , 0 . *50)- - 

j CALL SYMBOL (0. *-24.*H*H*F00T (26) *0. *50) 

XXa (XLENGTH-3,75) /2. 

CALL SYMBOL (XX*6T25*H*3lHWEIBULL DISTRIBUTION PARAMETERS. 0,0. 31 ) | 
999 RETURN ! 

!- ■ END 




APPENDIX E 


PDP PROGRAM TO CALCULATE ENDURANCE STRENGTH 
PARAMETERS FROM STAIRCASE TESTS 


C-FOCAL, 1969 

01.10 A "MINIMUM STRESS LEVEL", YP, ! 

01.20 A "STRESS INCREMENT", DP, ! 

01.30 A "NO OF SPECIMENS", NS, ! 

01.40 T "IF TEST IS BASED ON FAILURES THE CODE IS 0", ! 

01.50 T "IF BASED ON SUCCESSES TOE CODE IS 1", ! 

01.60 A "WHAT IS THE CODE?", Co, .' 

01.70 A "NO OF STRESS LEVELS IN TEST", I, !, ! 

01.74 S CU - 0 

01.75 S A * 0 

01.76 S B * 0 

01.77 T "NO OF SPEC IN EACH LEVL STARTING FROM TOE 2ND LOWEST", ! 

01.80 FOR J = 1, 1, 1-1; DO 4.0 

02.10 S SD = 1.62*DP*((NS*B-A+2)/NS 2+0.029) 

02.20 IF (CO) 2. 3, 2. 3, 2. 4 

02.30 S MU = YP+DP*(A/NS-.5) 

02.35 GOTO 2.7 

02.40 S MU = YP+DP*(A/NS- .5) 

02.70 T %10.03 "MEAN", MU, !, "STD DEV", SD, ! 

02.80 Q 

04.10 A NI, .' 

04.20 S CU = CU+1 

04.30 S A a A+CU*NI 

04.40 S B = B+CU 2*NI 




APPENDIX F 


PDP PROGRAM TO CALCULATE PAN WEIGHTS FOR DESIRED BENDING 179 

STRESS LEVELS FOR THE ANN ARBOR RESEARCH MACHINE 


C- FOCAL, 1969 

01.10 A "DIAMETER, D, 

01.11 A "LOWER STRESS", A, ! 

01.12 A "UPPER STRESS", B, ! 

01.13 A "INCREMENT", I, ! 

01.15 T "STRESS LEVEL"," MOMENT"," PAN LOAD", l 

01.20 F ST=A, I, B,; DO 2.0 

02.20 S M=3 . 1 416 *D+ 3*ST/ 32 

02.30 S P=((M+0.267)/4.07)-8.625 
02.40 T %8.40," " ,ST," ",M," 

0.350 GOTO 1.10 


m p t 
> r > • 
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